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SECTION  I 


INTRODUCTION 


1 . BACKGROUND 

It  is  a well  established  fact  that  the  structural  response  of 
hardened  weapon  system  facilities  to  nuclear  airblast  and  ground 
shock  loading  is  significantly  affected  by  the  manner  in  which 
the  structure  is  founded  upon  or  buried  within  the  soil  media. 
Although  there  is  an  acknowledged  awareness  that  the  structure 
interacts  with  the  adjacent  medium,  the  exact  nature  of  this  in- 
teraction is  imprecisely  defined. 

These  phenomena  and  related  aspects  have  generally  been  re- 
ferred to  as  soil  structure  interaction  (SSI)  . SSI  has  been  in- 
vestigated extensively  for  a number  of  years  from  different  points 
of  view,  and  the  determination  of  loads  on  buried  structures  is 
an  integral  part  of  the  subject. 

The  subject  of  SSI  constantly  recurs  in  connection  with  facil- 
ity response  problems  simply  because  there  is  no  logical  and  con- 
sistent way  of  avoiding  it.  The  dynamic  characteristics  of  a 
structure  in  contact  with  an  extended  material  (i.e.,  the  earth) 
differ  substantially  from  its  dynamic  characteristics  in  vacuo. 

The  surrounding  medium  has  the  effect  of  providing: 

• Resistance,  stiffening  the  structure  in  both  rigid  body 
and  flexible  modes  of  motion; 

• An  apparent  damping  of  the  structural  motions  by  trans- 
porting energy  in  the  form  of  a dynamic  disturbance  away 
from  the  structure  (i.e.,  radiation  damping); 

• Additional,  real  damping  due  to  the  dissipative  effects 
in  the  medium; 

• Fhock  mitigation  or  amplification  at  the  soil-structure 
interface  depending  on  the  properties  of  the  medium;  and 

• Arching  action  above  and  around  the  structure. 

In  terms  of  the  free  field,  the  inclusion  response  consists 
of  a perturbation  of  free  field  stress  waves  to  the  extent  that 
the  inclusion  displays  stiffness  or  mass  properties  which  differ 
from  the  material  for  which  the  free  field  stress  waves  are 


defined.  This  perturbation  also  displays  the  described  charac- 
teristic effects. 

Despite  previous  research  efforts  in  SSI,  it  remains  diffi- 
cult to  distinguish  and  identify  such  effects  in  real  structures. 
One  reason  for  this  is  that  in  most  real  structures  the  gross  mo- 
tions (rigid  body)  tend  to  follow  the  associated  free  field  dis- 
placement fields  with  relatively  small  displacement  perturbations. 
However,  the  fact  that  interaction  effects  are  relatively  small 
in  terms  of  displacement  fields  means  only  that  they  have  their 
greatest  influence  in  terms  of  velocities,  accelerations  and, 
structural  deformation  (i.e.,  stresses  and  strains)  which  are  more 
difficult  to  analyze  and  observe.  In  this  connection  note  that: 

• If  a structure  closely  follows  the  free  field,  this 
implies  that  there  is  enough  damping,  real  or  appar- 
ent, in  the  interaction  process  to  suppress  any  major 
oscillation  about  the  free  field  motions. 

• If  there  were  a strictly  imposed  condition  that  the 
structure  follow  the  free  field,  the  stress-strain  field 
throughout  the  whole  or  part  of  the  structure  would 

be  known  from  compatibility  and,  except  for  very  flex- 
ible structures,  these  would  be  completely  unrealistic. 

Although  the  subject  of  SSI  and  loads  on  buried  structures 
has  been  studied  to  some  extent  in  connection  with  such  applica- 
tions as  culvert  design  and  mining  operations,  early  attempts 
to  design  hardened  ground  facilities  gave  great  impetus  to  the 
field;  and  for  the  last  20  years  it  has  been  the  subject  of  ex- 
tensive research  in  connection  with  military  facilities  and, 
more  recently,  with  the  seismic  design  of  nuclear  power  plants. 

The  literature  associated  with  the  field  is  so  substantial  and 
widely  scattered  that  only  some  brief  comments  on  th'i  development 
of  the  field  with  an  emphasis  on  the  generation  of  sources  of 
data  and  the  application  of  such  data  to  the  design  process  are 
given  here.  References  1 through  57  are  cited  for  potential 
application  to  the  subject  investigation.  The  results  presented 
in  these  references  have  either  been  used  directly  or  reviewed 
for  possible  use  in  this  research  program.  Several  attempts  at 
a detailed  literature  review  have  been  made  including  those  in 
references  1 through  6. 
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Numerous  research  investigations  have  been  undertaken  to 
study  and  quantify  various  phenomena  associated  with  buried 
structures  under  dynamic  loads  in  connection  with  the  need  to  de- 
sign or  harden  military  facilities  against  nuclear  attack.  These 
early  studies  employed  data  from  nuclear  tests,  small-scale  lab- 
oratory experiments,  existing  solutions  for  classical  inclusion 
problems,  and  a design  rationale  largely  adopted  from  soil  me- 
chanics. Early  centers  for  studies  of  this  type  were  at  IIT  Re- 
search Institute  (IITRI) , then  Armour  Research  Foundation,  the 
University  of  Illinois,  and  the  Massachusetts  Institute  of  Tech- 
nology. This  phase  of  research  culminated  with  the  appearance 
of  a variety  of  handbooks  and  technical  summaries  such  as  the 
early  Air  Force  Design  Manual  (Ref.  5),  the  Compendium  of  Air- 
blast  Effects  (Ref.  6)  and  a series  of  Corps  of  Engineers  manuals 
entitled  Design  of  Structures  to  Resist  the  Effects  of  Atomic 
Weapons  (Ref.  7).  These  documents  generally,  and  the  Air  Force 
Design  Manual  in  particular,  took  a conventional  and  somewhat 
conservative  approach  to  the  question  of  loads  on  buried  struc- 
tures. The  design  process  employed  was  basically  quasi-static 
and  the  structure  was  essentially  designed  for  the  estimated  free 
field  pressure  wave.  The  beneficial  effects  of  soil  arching  above 
and  around  the  structure  were  also  taken  into  account;  however, 
the  possibility  of  load  magnification  on  the  structure  due  to 
reflection  of  the  soil  stress  wave  at  the  interface  was  generally 
dismissed  on  the  grounds  that  little  or  no  evidence  for  such  phe- 
nomena could  be  found  in  the  available  test  data. 

Subsequent  activity  saw  the  decreased  reliance  on  full-scale 
test  data  (the  moratorium  on  aboveground  testing  had  come  into 
effect),  the  increased  use  of  medium-scale  field  simulations,  and 
especially,  the  increased  use  of  large-scale  computations  and  com- 
puter programs . Early  computer  work  involved  the  extension  of 
more  or  less  classical  approaches  (i.e.,  series  and  modal  solu- 
tions) as  in  references  8,9  and  10,  and  the  development  of  so-called 
SSI  models  in  which  effects  at  the  structure-soil  interface  were 
represented  by  stiffness  and  damping  terms.  An  early  model  of  this 
type  was  developed  at  IITRI  (Ref.  11,  12,  and  13)  in  which  the 


normal  stress  ag  (i.e.,  unit  load)  on  the  surface  of  the  structure 
was  represented  as 

°s  = °f+k(uf-us)+s(uf-us)  (1) 

where  a ^ is  the  normal  stress  in  the  free  field  (i.e.,  in  the  soil 
if  the  structure  were  not  present) ; U£  and  ug  are  corresponding 
normal  displacements  of  the  free  field  and  the  structure;  and  u^ 
and  ug  are  the  velocities.  Such  representations  were  for  a time 
extensively  used  by  IITRI  and  others  in  the  evaluation  and  design 
of  facilities  including  the  MINUTEMAN  and  SAFEGUARD  system  struc- 
tures (see  Ref.  14  through  19  for  instance).  The  approach  was 
developed  as  a generalization  of  fairly  old  techniques  related  to 
culvert  and  tunnel  design;  and  while  it  has  fallen  into  disuse  in 
recent  years  due  to  theoretical  objections  and  the  development  of 
finite  element  and  other  techniques,  it  remains  a potentially  val- 
uable tool  for  the  evaluation  of  loads  on  buried  structures  and 
known  classical  solutions  can  be  set  up  in  this  form.  The  method 
displays  several  known  ingredients  of  the  inclusion  response  solu- 
tion, namely,  the  explicit  dependence  on  the  free  field  stress, 
the  recognition  that  the  surrounding  medium  has  stiffness,  and  the 
incorporation  of  wave  dispersion  or  interface  reflection  in  the 
damping  term. 

The  solution  to  classical  inclusion  problems  has  been  the  sub- 
ject of  numerous  investigations  throughout  the  most  active  period 
of  SSI  research.  The  response  to  spherical,  elastic  inclusions  to 
steady  state  dilatational  wave  strains  has  been  studied  with  re- 
spect to  stress  concentrations  in  the  surrounding  medium  (Ref.  20) 
and  scattering  phenomena  (Ref.  21).  In  both  cases  the  numerical 
results  are  few  and  apply  principally  to  the  limiting  cases  of  rig- 
id and  hollow  spheres.  The  analogous  acoustic  inclusion  problem 
(i.e.,  acoustic  sphere  in  acoustic  medium)  has  been  solved  in  refer- 
ence 22.  The  rigid  spherical  inclusion  has  also  been  studied  by 
C.  C.  Mow  in  references  23  through  25,  with  emphasis  on  the  inclu- 
sion response  for  a transient  compressive  pulse  with  exponential 
decay.  These  papers,  especially  reference  24,  contain  the  principal 
qualitative  information  derivable  from  the  classic  inclusion  lit- 
erature in  terms  useful  for  the  structural  inclusion  problem. 
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Briefly  these  findings  are: 
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• The  velocity  of  the  inclusion  exceeds  the  incident 
free  field  particle  velocity  only  when  the  inclusion 
is  heavier  than  the  same  volume  of  the  surrounding 
medium. 

• The  solution  contains  a free  vibration  term  which  de- 
cays exponentially  as  energy  is  carried  away  from  the 
inclusion  response  by  wave  action. 

• The  velocity  and  acceleration  response  of  the  inclu- 
sion may  differ  substantially  from  the  corresponding 
free  field  motions  even  though  the  differences  in 
displacement  components  are  not  great. 

The  first  of  these  observations  is  often  cited  as  an  argument 
that  there  can  be  no  amplification  in  the  gross  motion  of  a struc- 
tural inclusion  since  these  are  inherently  hollow,  lined  cavi- 
ties and  therefore  of  lower  average  density  than  the  surrounding 
media . 

The  corresponding  two-dimensional  inclusion  problems  have 
been  more  extensively  studied  although  t,hey  generally  remain  lim- 
ited to  elastic  materials  and  simple  geometries.  References  26 
through  31  are  concerned  with  the  impingement  of  plane  waves  on 
cylindrical  cavities.  Stress  concentrations  are  obtained  in  refer 
ence  30  for  harmonic  dilatational  wave  trains  as  functions  of 
wavelength  and  Poisson's  ratio.  References  27,  28.  and  29  consider 
the  response  of  cylindrical  cavities  for  transient  compression 
pulses  and  generate  influence  functions  which  are  subsequently 
used  in  obtaining  the  response  of  a cylindrical  shell  used  as 
a cavity  liner  (Ref.  32  and  33).  The  equivalent  problem  for  a 
cylindrical  shell  in  an  acoustic  medium  is  given  in  references 
34  and  35. 

The  rigid  circular  inclusion  is  treated  for  harmonic  dilata- 
tional waves  in  references  36  and  37,  with  reference  36  providing 
information  on  the  relation  of  inclusion  response  to  wavelength 
and  inclusion  mass  (similar  to  that  given  in  Ref.  24) , and  refer- 
ence 37  giving  stress  concentration  data.  Exact  elasticity  solu- 
tions for  thick  cylindrical  inclusions  are  given  in  reference  38 
for  harmonic  dilatational  waves  and  in  reference  39  for  incident 
shear  waves.  References  40  and  41  consider  the  response  of 
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cylindrical  cavities  and  shells  to  plane  waves  which  are  obliquely 
incident  (i.e.,  propagate  along,  in  addition  to  across,  the  axis 
of  the  cylinder) . 

A class  of  problems  that  can  be  called  the  surface  inclusion 
or  footing  problem  involves  the  response  of  various  types  of  slabs 
resting  on  the  surface  of  a semi-inf inite , and  perhaps  layered, 
elastic  medium.  These  differ  from  the  spherical  and  cylindrical 
classical  problems  only  in  that  the  slab  does  not  replace  any  of 
the  surrounding  medium.  The  classic  papers  for  this  problem  are 
by  Bycroft  (Ref.  42  and  43)  and  deal  with  the  forced  harmonic  vi- 
bration of  a rigid,  circular  surface  footing  in  four  distinct 
modes.  The  solutions  consist  of  driving  force  versus  frequency, 
amplitude,  and  dimensionless  constants  and  have  proved  useful  in 
the  analysis  of  the  response  of  surface  structures.  Reference  44 
provides  similar  solutions  for  the  rectangular  slab  and  indicates 
that  these  results  are  similar  over  a range  of  footing  aspect 
ratios  to  those  for  the  circular  slab.  Unlike  fully  buried  struc- 
tures, heavy  surface  structures  can  display  significant  amplifi- 
cation effects  since  they  are  not  limited  by  a displaced  volume 
of  soil  and  can  generate  significant  surface  mass  densities. 

More  recent  activity  in  the  SSI  research  has  seen  the  increased 
reliance  on  large-scale  computations  and  simulations  of  the  re- 
sponse of  buried  structures.  This  is  due  almost  entirely  to  the 
emergence  of  the  finite  element  analysis  technique  as  a powerful 
tool  making  such  simulation  possible.  Although  a wide  variety  of 
programs  and  techniques  have  been  used  for  the  analysis  of  buried 
structures,  some  programs  which  have  been  developed  and  used  spe- 
cifically for  this  purpose  include  INDEPS  (Ref.  45),  NOFEAR  (Ref.  46) 
and  DYNAX  (Ref.  47),  and  SAMSON  (Ref.  48);  the  last  two  were  de- 
veloped at  IITRI . 

It  is  noteworthy  that  throughout  these  later  developments  com- 
paratively little  effort  has  been  devoted  to  extract  design  type 
data  on  loads  and  response  or  to  the  preparation  of  summaries  and 
systematic  simplified  procedures.  The  new  edition  of  the  Air 
Force  Design  Manual  (Ref.  49)  presents  little  new  data  relating 
to  the  loads  on  buried  structures  over  that  which  were  available 
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in  the  earlier  editions.  In  fact,  the  i manual  treats  tnis 
phase  of  design  almost  exclusively  in  terms  of  response  of  rigid 
structures  with  soil  response  represented  by  dissipative  terms  of 
the  form  pcV  corresponding  to  plane  elastic  wave  propagation.  This 
is  equivalent  to  the  simplest  form  of  the  SSI  model  described  in 
equation  (1)  with  k = 0 and  s=pc. 

2.  OBJECTIVES 

This  research  project  has  been  directed  toward  the  develop- 
ment of  simplified  techniques  for  determining  loads  on  buried 
structures  of  cylindrical  geometry  subjected  to  severe  dynamic 
blast  load  environments.  The  data  forming  the  basis  for  these 
design  loads  have  been  primarily  derived  from  numerical,  finite 
element  analyses  for  a wide  range  of  loading,  soil,  and  structural 
parameters.  The  numerical  results  presented  herein  have  been 
generated  exclusively  from  this  research  investigation. 

The  primary  objective  of  this  program  has  been  to  generate 
SSI  loads  and  internal  cylinder  forces  for  horizontally  oriented 
cylinders.  Both  static  and  dynamic  solutions  have  been  obtained 
in  order  to  study  in  a more  or  less  systematic  manner  the  SSI 
loads  and  internal  cylinder  forces  for  variations  of  the  struc- 
tural stiffness,  depth  of  burial,  and  soil  media  parameters. 

During  the  course  of  this  project,  a concept  was  postulated 
that  air  blast  induced  SSI  loads  could  be  represented  by  only  a 
few  components  of  the  modal  representation  of  the  SSI  loads.  For 
example,  a buried  cylinder  in  an  infinite  elastic  medium  is  sub- 
jected to  only  the  uniform  (mode  0)  and  ovaling  (mode  2)  load 
components  for  a statically  applied  biaxial  state  of  free  field 
stress.  With  this  concept  in  mind,  various  static  and  dynamic 
buried  structure  analyses  were  performed  to  determine  the  appro- 
priateness of  characterizing  SSI  loads  by  a few  modes  for  a wide 
range  of  structure,  soil,  depth  of  burial  and  loading  parameters. 

3.  SCOPE  OF  WORK 

Attempts  have  been  made  throughout  the  project  to  cover  a 
rather  wide  range  of  loading,  structure,  and  soil  media  parameters. 
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However,  the  depths  of  burial  of  the  structure  that  were  consid- 
ered are  relatively  shallow  since  it  has  been  postulated  that 
ground  surface  effects  can  be  generally  ignored  for  depths  of 
burial  greater  than  two  times  the  cylinder  diameter.  Most  of  the 
results  are  presented  in  nondimensional  form;  however,  if  any  spe- 
cific dimensions  were  required,  such  as  the  diameter  of  the  cylin- 
der, then  a 4 meter  inside  diameter  has  been  employed.  The  cylin- 
der has  also  been  idealized  throughout  the  study  as  an  elastic 
material  with  thicknesses  and  material  parameters  corresponding  to 
concrete . 

The  initial  phase  of  the  project  was  devoted  to  a brief  review 
of  the  literature  which  was  summarized  in  section  1.1.  Other  ac- 
tivities required  early  in  the  project  were  to  establish  a commu- 
nication link  with  the  AFWL  computer  facilities  at  Kirtland  AFB, 

New  Mexico.  This  effort  was  required  since  the  majority  of  the 
numerical  calculations  were  to  be  performed  on  the  CDC  6600  and 
7600  computers  at  Kirtland  AFB  through  a remote  batch  configura- 
tion via  a dial-up  terminal  at  IITRI  in  Chicago. 

Several  simplified  methods  of  analysis  for  determining  loads 
on  buried  structures  have  been  developed  in  the  past.  Cne  such 
technique  that  is  based  on  experimental  results,  analytical  inves- 
tigations, and  engineering  judgement  was  defined  early  in  this  proj- 
ect by  N.  M.  Newmark  Consulting  Engineering  Services  and  is  pre- 
sented in  appendix  II.  This  method  does  not  reflect  the  numerical 
results  of  the  analyses  conducted  in  this  investigation. 

The  early  numerical  results  were  devoted  to  the  generation  of 
some  elastic  material,  static  load  solutions  for  variations  of  the 
cylinder  and  soil  media  parameters.  The  results  of  these  studies 
are  given  in  section  II. 

There  are  some  indications  that  the  principal  components  of 
SSI  loading  on  a horizontally  oriented  cylinder  can  be  represented 
by  a relatively  few  modal  components  of  the  loading.  This  simpli- 
fied loading  concept  appeared  to  be  plausible  since  the  uniform 
component  (mode  0)  and  ovaling  component  (mode  2)  are  the  only 
sources  of  loading  for  an  elastic  cylinder  embedded  in  an  isotropic 
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elastic  media  of  infinite  extent  subjected  to  a free  field  biaxial 
state  of  static  loading.  The  major  effort  undertaken  during  the 
project  was  devoted  to  the  generation  of  dynamic  finite  element 
numerical  results  in  an  effort  to  ascertain  whether  or  not  a few 
modes  can  be  selected  which  will  adequately  represent  the  loading 
on  a structure.  The  SAMSON  computer  code  was  used  for  these  anal- 
yses to  obtain  solutions  to  a basic  set  of  structure,  soil,  and 
loading  parameters.  A total  of  32  cases  was  defined  to  represent 
what  was  considered  to  be  reasonable  extreme  values  of  the  param- 
eters that  were  varied.  The  five  parameters  considered  were: 

• Radius  to  thickness  ratio  of  cylinder 

• Depth  of  burial 

• Peak  airblast  overpressure 

• Duration  of  overpressure  loading 

• Soil  types  (nonlinear  material) 

Two  additional  elementary  cases  were  also  solved  using  nonlinear 
reinforced  concrete  cylinder  parameters.  Also,  two  solutions  were 
obtained  using  an  elastic  soil  material  in  conjunction  with  the 
elastic  concrete  tunnel  liner.  These  finite  element  SAMSON  solu- 
tions and  the  manner  in  which  the  numerical  results  are  processed 
and  displayed  are  presented  in  section  III. 

Section  III  also  contains  the  method  of  analysis  that  was  used 
to  judge  the  postulation  that  only  a few  of  the  modal  components 
of  the  load  are  required  to  adequately  represent  the  essential 
SSI  loads.  A discussion  of  the  results  of  the  SAMSON  analyses  and 
the  problems  encountered  in  determining  which  modal  load  components 
are  important  is  given  in  section  III.  Conclusions  and  recommen- 
dations for  future  work  are  in  section  IV. 

Appendix  A gives  the  results  of  a number  of  free  field  solu- 
tions that  were  generated  involving  the  two  soils  considered  and 
the  four  airblast  loading  waveforms.  Appendix  B presents  a simple 
method  for  determining  loads  on  buried  structures  that  was  defined 
by  N.  M.  Newmark  Consulting  Engineering  Services  at  the  beginning 


I « 

r 

of  this  project.  Appendix  C is  presented  to  indicate  the  mag- 
nitude of  errors  that  can  be  expected  if  the  stresses  in  the  soil 
elements  surrounding  the  structure  are  used  to  determine  the  SSI 
loads . 

Volume  II  of  this  report  contains  the  bulk  of  the  SAMSON  solu- 
tions for  each  of  the  cases  that  were  run.  These  results  include: 
(1)  tables  of  maximum  and  minimum  values  of  the  cylinder  internal 
forces  and  SSI  stresses  around  the  liner,  (2)  modal  amplitude  his- 
tories of  the  cylinder  forces  and  the  SSI  stress,  (3)  time  history 
plots  of  the  cylinder  moment  and  thrust  and  the  SSI  radial  stress 
at  the  crown,  springline,  and  invert  points. 


SECTION  II 
STATIC  ANALYSES 

This  section  presents  the  results  of  several  static-elastic 
analyses  that  were  performed  to  investigate  the  effects  of  depth 
of  burial  and  structure  and  soil  media  stiffness  on  the  SSI  stresses 
and  the  cylinder  internal  forces.  These  investigations  are  limi- 
ted to  plane  strain  analyses  of  horizontally  oriented  cylinders, 
and  hence  they  are  applicable  to  the  analysis  of  a long  cylinder 
where  the  effects  of  the  ends  of  the  cylinder  can  be  neglected. 

Both  the  structure  and  the  soil  media  have  been  indealized  as  linear 
elastic  materials. 

1.  INFINITE  MEDIA  ANALYSES 

These  analyses  consider  an  elastic  liner  embedded  in  a homo- 
genous elastic  soil  media  of  infinite  extent.  A biaxial  free  field 
state  of  stress  in  the  soil  is  assumed,  and  the  interface  between 
the  soil  and  structure  can  be  idealized  as  either  free  to  slip  (i.e., 
no  interface  friction)  or  fully  bonded  (i.e.,  infinite  interface 
friction) . This  problem  has  been  treated  extensively  in  the  past 
by  Burns  and  Richards  (Ref.  50)  and  Hoeg  (Ref.  51  and  52) . Burns 
and  Richards  restricted  their  investigation  to  the  case  where  the 
biaxial  free  field  stress  state  is  derived  from  a uniaxial  state 
of  strain,  and  solutions  are  given  for  the  calculation  of  the  media 
stresses  and  the  cylinder  force  resultants.  The  major  disadvantage 
of  the  uniaxial  strain  case  is  that  the  biaxial  stress  ratio  is 
dependent  on  Poisson's  ratio  of  the  soil. 

Hoeg's  studies  are  not  restricted  to  a uniaxial  strain  case, 
and  the  biaxiality  ratio  can  be  arbitrarily  specified.  However, 
these  investigations  do  not  include  solutions  for  internal  force 
resultants  of  the  cylinder  and  a typographic  sign  error  exists  in 
one  of  the  equations  which  can  lead  an  unsuspecting  user  to 
erroneous  conclusions. 

* The  terms  cylinder,  tunnel,  and  liner  are  used  interchangeably 
to  denote  the  right  circular  cylinder  buried  structure. 
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The  following  equations  are  based  on  appropriately  modified  and, 
where  necessary,  corrected  results  from  reference  51. 

2.  SOIL  MEDIA  STRESS  ANALYSIS 

Referring  to  figure  1,  the  nondimens ional  media  stresses  at 
any  nondimens ional  polar  coordinate  in  the  media  (p,9)  are  given 
by 


Radial  Stress,  ar  = 

ar/p  = arQ  + or2  cos  29 

(2) 

Circumferential  Stress, 

°0  = ag/p  = aQO+^02  cos  20 

(3) 

Shear  Stress,  Trg 

= Tr97P  = Tr92  sin  20 

(4) 

where 


k = biaxiality  ratio  of  free  field  soil  stress  and  the  sub- 
scripts 0 and  2 refer  to  the  uniform  and  ovaling  components. 


p = R/r  (<  1.0)  (5a) 

orQ  = -(l+k)(l-alP2)/2  (5b) 

ar2  = -(l-k)(l-3a2  p4-4a3  p2)/2  (5c) 

0go  = -(l+k)(l+a1  p2)/2  (5d) 

aQ2  = (l-k)(l-3a2  p4)/2  (5e) 

Tr02  = (l-k)(l+3a2  p4+2a3  p2)/2  (5f) 


The  coefficients  a-^,  a2,  a3  have  been  determined  for  two  ideal- 
ized soil-structure  interface  conditions;  namely,  the  fully  bonded 
or  no  slip  (NS)  case  and  the  zero  friction  or  free  slip  (FS)  case. 
For  the  NS  case  the  coefficients  are 


a = (l-2v)(C-l) 
1 (1-2 v)  c+1 


i2  = (l-2v)(l-C)  F-(l-2v)2  C/2+2  /b 


-{[ 


l+(l-2v)  c 


F - (l-2v)  C/2-2  y 


/b 


where  b = [ (3-2v)  + (i-2  v)  c]f  + (2 . 5-8v  + 6v2)  C + 6-8' 


(ba) 

(6b) 

(6c) 

(7a) 


and 


C = a(R/t)/ (1- v) , compressibility  ratio 
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(7b) 


F = 2C  (l-2v) (R/ 1)  , flexibility  ratio  (7c) 

★ 

a = M /Ec , material  modulus  ratio  (7d) 

M*  = w-V^-n-  - , , constrained  modulus  of  soil  (7e) 

(1+v) (I-ZvJ ’ 

E,v  = modulus  of  elasticity  and  Poisson's  ratio  of  soil  media 
_ 2 

Ec  = Ec/(l-vc),  apparent  modulus  of  cylinder  (7f) 

Ec,vc  - modulus  of  elasticity  and  Poisson's  ratio  of  cylinder 

For  the  FS  case  the  coefficients  are 

a,  = (1-"2v)  (8a) 

1 (l-2v)  C+l 


a2 


2F+1  - 2v 

2f+5  - 6v 


a3 


2F-1 


If  the  free  field  is  in  a state  of  uniaxial  strain,  i.e.,  only 
the  vertical  strain  is  nonzero,  then  the  biaxiality  ratio  is  de- 
termined from  Poisson's  ratio  of  the  soil  as 

k = v/(l-v),  uniaxial  strain  case  (9) 


Equations  (2)  through  (8)  for  the  calculation  of  soil  stress 
are  independent  of  equation  (9),  and  the  biaxiality  ratio  can  be 
any  value  desired. 

The  solutions  for  the  FS  and  the  NS  idealizations  have  been 
used  to  compute  the  radial  and  circumferential  stress  variation 
in  the  vicinity  of  the  cylinder  along  the  crown  ray  (0  = 0)  and 
the  springline  ray  (9  = 90  deg)  . Since  the  solution  is  for  the 
unbounded  or  infinite  media  case,  the  invert  ray  (9  = 180  deg) 
stress  variation  is  the  same  as  for  the  crown  ray.  Results  were 
obtained  for  several  cylinder  stiffnesses  and  a typical  soil 
media  with 


14 


E = 25  ksi  (172.4  MPa)  or  M*  = 30  ksi  (206.9  MPa) 

v = 0.25 
k = 1/3 

These  soil  and  structure  parameters  will  be  referred  to  as  Case  A. 
The  cylinder  material  parameters  are  appropriate  for  a concrete 
tunnel  liner,  viz., 

Ec  = 3000  ksi  (20,700  MPa) 

v = 0.2 

c 

Results  are  shown  in  figure  2 for  four  cylinder  geometries  rep- 
resented by  radius  to  thickness  ratios  (R/t)  of  15,  8,  6 and  4. 

For  a tunnel  with  a 6.6  ft  (2  m)  inside  radius,  these  R/t  ratios 
correspond  to  cylinder  thicknesses  ranging  from  approximately  0.5 
to  2 ft  (0.152  to  0.610  m) . For  this  figure,  a rather  unconven- 
tional method  has  been  used  to  display  the  radial  and  circumferen- 
tial soil  stress  variations  in  a condensed  manner.  Note  that  all 
of  the  stresses  are  compression  (negative)  except  the  circumferen- 
tial stress  near  the  crown  for  the  FS  case.  The  compressive  (-) 
radial  stresses  have  been  displayed  above  the  abscissa  (crown  or 
springline  ray);  whereas,  the  compressive  (-)  circumferential 
stresses  are  displayed  below  the  abscissa.  Hence,  the  tension  (+) 
circumferential  stresses  at  the  crown  are  plotted  above  the 
abscissa . 

An  examination  of  the  curves  in  this  figure  indicates  that  the 
soil  stresses  are  very  nearly  equal  to  the  free  field  stress  for 
points  approximately  1 to  1.5  diameters  away  from  the  edge  of  the 
liner.  These  curves  also  show  that  the  soil  stresses  adjacent  to 
the  cylinder  are  significantly  different  for  the  NS  and  the  FS 
solutions,  especially  the  circumferential  stresses.  The  NS  case 
is  considered  a more  realistic  occurrence  than  the  FS  case.  For 
example,  it  has  been  observed  that  when  highway  tunnel  liners  that 
have  been  in  service  for  some  time  are  excavated,  soil  particles 
are  attached  or  bonded  to  the  entire  surface  of  the  liner.  This 
indicates  that  if  relative  tangential  motion  occurs  between  the 
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structure  and  the  soil,  it  takes  place  in  the  soil  because  the 
shear  capacity  of  the  soil  is  exceeded  and  not  because  there  is 
low  friction  strength  at  the  soil-structure  interface. 

3.  CYLINDER  FORCE  RESULTANT  ANALYSIS 


The  internal  force  resultants  (moment  and  thrust)  of  the  cylin- 
der can  be  computed  in  the  following  manner.  The  elastic  law  for 
the  cylinder  moment  is  (Ref.  53,  p 211) 


M = 


K 

P 


(10a) 


where 

K = Ec  t3/12  (10b) 

The  radial  displacement  of  the  cylinder  is 


w = 0 . 5p 


_R 

M* 


(1+k) (1-v) 


1 + 


l-2v 


(1-k)  l + a2  + 4(l-v)  a3Jcos  2o| 


(10c) 


The  nondimens ional  radial  displacement  of  the  cylinder  can  be 
written  as 


w = wM*/ (pR)  = w0  + w2  cos  (10d) 

where 

wQ  = 0 . 5 (1+k)  (1-v)  (l  + T^J  (lOe) 

w2  = 0 . 5 (1-k)  j^l  + a2  + 4(l-v)a3j  (lOf) 

From  equations  (10a)  and  (10b)  the  nondimens ional  cylinder  moment 
is  obtained  as 

M = M/pR2  = Mq+M2  cos  20  (lOg) 

where 

M0  = -C1(l+k)  (1-v)  (l  +T~2v|  dOh) 
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M2  = 3C1(l-k)  l + a2  +4(l-v)a3 

C1  = = ^ ^ 


(10i) 


(10j) 


The  differential  equation  of  radial  equilibrium  for  the  cylinder 
is  (Ref.  53,  p 211) 


d*M  tjt  — n TJ 2 

— ir  - RT  = a R 
d9Z  r 


The  radial  stress  (ar)  referred  to  in  this  equation  is  the  inter- 
face SSI  stress,  computed  from  equations  (2)  and  (5)  with 
P = 1.  Performing  the  indicated  operations  in  equation  (11)  and 
solving  for  the  nondimens ional  thrust  yields 


T = T/pR  = Tq  + T2  cos  29 


where 


T0  “°r0 


@ P=1 


= (l+k)(l-ai)/2 


(Ha) 


(Hb) 


T2  = -°r2 


-4  M2 

J@  P=1 


= 0.5  (k-1)  jl-3a2-4a3-24  C1(]T^)  l+a2+4(l-v) a3  | 


(He) 


A very  close  approximation  to  the  nondimensional  cylinder 
force  modal  amplitudes  can  be  written  as  a function  of  the  non- 
dimensional  SSI  modal  stresses  as 


m _ ro 
M = 

° ' 12  (R/t) 


(lid) 


M2  = ("0r2  + Tr02/2)/3 


(He) 


To  " aro 


(Ilf) 


T2  (ar2  " 2 Tr02)/3 


(llg) 


l 


The  following  four  fundamental  nondimensional  parameters 
are  sufficient  to  completely  specify  the  material  and  geometric 
characteristics  of  the  problem: 

a,  modulus  ratio,  equation  (7d) 

k,  biaxiality  ratio  of  free  field  stress 

v,  soil  Poisson's  ratio 

R/t,  cylinder  radius  to  thickness  ratio 

In  an  effort  to  study  the  nondimensional  variation  of  the  SSI 
stresses  (®r>^rg)  an<^  the  cylinder  internal  forces  (M,T)  around 
the  cylinder  (9  = 0 to  90  deg)  , a number  of  elementary  solutions 
have  been  computed  using  the  appropriate  equations  (2)  through 
(11)  for  the  NS  case  at  the  interface  idealization. 

As  before,  a typical  soil  medium  postulated  (E  = 25  ksi  = 

172  MPa,  v = 0.25,  k = 1/3)  and  a concrete  cylinder  is  assumed 
(Ec  = 3000  ksi  = 20,690  MPa,  = 0.2)  with  the  same  four  radius- 
to- thickness  ratios.  For  this  case,  the  fundamental  parameters  are 

Case  A:  a = 0.0096 
k = 0.333 
v = 0.25 

R/t  = 15,  8,  6,  4 

The  circumferential  variations  of  the  nondimensional  SSI 
stresses  and  the  cylinder  internal  forces  for  this  case  are  shown 
in  figure  3. 

Solutions  to  three  additional  cases  have  also  been  obtained 
to  indicate  the  influence  of  the  parameter  a on  the  results.  For 
these  cases,  the  biaxiality  ratio  k has  been  taken  as  zero  so  that 
the  results  can  be  used  to  obtain  solutions  for  any  biaxiality 
ratio.  This  can  be  done  by  transforming  the  k = 0 results  by 
90  deg,  multiplying  by  the  desired  biaxiality  ratio,  and  adding 
this  solution  to  the  original  k = 0 solution.  Only  the  param- 
eter a has  been  varied  between  these  three  cases;  the  remaining 
fundamental  parameters  for  each  case  have  been  taken  as 
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k = 0 
v = 0.25 

R/t  = 15,  8,  6,  4 


1 

1 


cases  B,  C and  D 


The  modulus  ratio  and  figures  containing  the  SSI  stresses  and 
cylinder  internal  forces  for  each  case  are: 

Case  BCD 

a 0.0096  0.096  0.96 

figure  456 

Case  B has  the  same  value  of  a as  case  A,  but  assumes  there  is  no 
horizontal  free  field  stress.  Cases  C and  D use  values  of  a that 
are  one  and  two  orders  of  magnitude  greater  than  case  B.  Thus,  if 
the  tunnel  liner  material  is  again  assumed  to  be  concrete,  then  the 
constrained  modulus  of  the  soil  media,  M*,  would  be  30,  300  and 
3000  ksi  (206.9,  2069  and  20,690  MPa)  for  cases  B,  C and  D,  re- 
spectively. Case  D is  representative  of  a very  competent  rock- 
like soil  media,  and  case  C would  represent  a material  intermediate 
between  a typical  soil  and  rock-like  media. 

Table  I gives  the  mode  0 and  mode  2 component  amplitudes  of 
the  nondimens ional  SSI  stresses  and  the  cylinder  internal  forces 
for  each  of  the  four  sets  of  material  parameters  (cases  A,  B,  C 
and  D) . These  amplitudes  have  been  computed  for  each  of  the 
structure-soil  media  interface  idealizations,  i.e.,  the  fully 
bonded  or  NS  assumption  and  the  zero  friction  or  FS  condition. 

The  results  of  cases  B,  C and  D demonstrate  a number  of  in- 
teresting things  such  as  the  arching  action  through  the  soil 
around  the  liner  for  the  more  flexible  cylinders  of  case  B and  very 
high  magnitudes  of  arching  around  the  liner  for  all  cylinders  con- 
sidered for  cases  C and  D.  Also,  the  SSI  shear  stress  (xrg)  and 
the  cylinder  thrust  (T)  are  not  terribly  sensitive  to  the  cylinder 
stiffness  (R/t)  or  the  material  modulus  ratio  (a)  until  the  soil 
media  modulus  approaches  that  of  the  _ylinder.  As  expected,  the 
bending  moments  in  the  liner  (M)  are  very  sensitive  to  both  the 
modulus  ratio  (a)  and  the  cylinder  stiffness  (R/t). 
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Table  I 

NONDIMENS IONAL  SSI  AND  CYLINDER  FORCE  AND  DISPLACEMENT  COMPONENTS 
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Note:  The  mode  0 component  for  the  FS  interface  condition  Is  the  same  value  as  for  the  NS  condition. 


4.  DEPTH  OF  BURIAL  ANALYSIS 

This  portion  of  the  static  analysis  investigations  again  con- 
siders an  elastic,  horizontally  oriented  cylinder  embedded  in  a 
homogeneous  elastic  soil  media.  The  influence  of  depth  of  burial 
of  the  cylinder  was  studied  by  performing  static  finite  element 
analyses  for  ratios  of  depth  of  burial  to  cylinder  radius  (h/R) 
values  of  0.5,  1.0,  2.0  and  4.0.  The  SAP4  computer  code  (Ref.  54) 
was  employed  for  all  of  these  investigations.  The  following  soil 
media  and  structural  parameters  were  used  for  each  of  the  four 
depth  of  burial  finite  element  investigations: 

Soil  Parameters:  E = 25  ksi  (172.4  MPa) 

v = 0.25 

Cylinder  Parameters:  E = 3000  ksi  (20,700  MPa) 

v =0.2 
c 

R/t  = 15,  8,  6,  4 

The  ground  surface  is  subject  to  a uniform  pressure  loading 
(p)  as  indicated  in  figure  1.  Since  the  postulated  problem  is 
symmetric  about  a vertical  plane  through  the  axis  of  the  cylinder, 
only  half  of  the  system  needs  to  be  modeled.  A plane  strain  anal- 
ysis has  been  employed  for  the  soil  media  for  each  of  the  SAP4 
numerical  solutions.  The  results  of  the  infinite  media  analyses 
illustrated  in  figure  2 imply  that  for  a static  analysis  the  in- 
fluence of  the  cylindrical  inclusion  is  essentially  negligible 
for  points  greater  than  1.5  times  the  cylinder  diameter.  This 
conclusion  is  based  on  the  observation  that  the  state  of  stress 
for  points  1.5  diameters  from  the  cylinder  is  essentially  the  same 
as  the  free  field  state  of  stress.  Therefore,  the  finite  element 
mesh  boundaries  (vertical  "lateral"  and  horizontal  "lower")  were 
located  at  two  diameters  from  the  edge  of  the  cylinder.  The 
boundary  conditions  imposed  on  the  vertical  and  lower  horizontal 
finite  element  mesh  perimeters  were  such  that  the  nodes  were  free 
to  translate  parallel  to  the  boundary  but  were  restrained  from 
all  motion  normal  to  the  boundary  and  were  also  not  permitted  to 
rotate.  These  conditions  simulate  a roller  bearing  surface  at  the 
boundaries.  The  remaining  grid  points,  including  the  upper  ground 
surface  nodes,  were  unrestrained. 


A considerable  effort  was  made  to  keep  the  aspect  ratio  of 
each  element  close  to  1 for  each  of  the  finite  element  meshes 
used.  The  resulting  grid  works  are  shown  in  figures  7 through  10 
for  the  four  depth  of  burial  h/R  values  of  0.5,  1.0,  2.0  and  4.0, 
respectively.  As  previously  indicated,  two-dimension  plane  strain 
elements  were  used  to  model  the  soil  media  (type  4 elements  - 
SAP 4) . 

Twenty-five  equally  spaced  node  points  were  used  to  simulate 
the  half-model  of  the  soil-structure  interface.  Thus,  each  pair 
of  adjacent  nodes  is  located  7.5  deg  apart  at  the  soil-structure 
interface.  The  tunnel  liner  was  modeled  by  interconnecting  24 
beam  "line"  elements  (type  2 elements  - SAP4)  to  these  25  node 
points.  Therefore,  the  finite  element  mesh  simulates  the  NS  condi- 
tion at  the  structure  media  interface. 

Each  of  the  four  depth -of -burial  meshes  (figures  7 through  10) 
was  used  to  obtain  numerical  solutions  for  tne  four  cylinder 
stiffnesses  (R/ t = 15,  8,  6,  4);  thus,  solutions  were  obtained  for 
a total  of  16  cases.  The  cylinder  stiffness  (R/t)  was  simulated 
by  specifying  appropriate  moments  of  inertia  and  cross-sectional 
areas  for  the  beam  elements.  The  number  of  nodes  and  media  ele- 
ments used  for  each  of  the  finite  element  meshes  is  indicated 
below. 

h/R  Nodes  Elements 


0.5 

173 

159 

1.0 

186 

177 

2.0 

190 

188 

4.0 

202 

198 

The  results  from  the  16  SAP4  numerical  solutions  are  summar- 
ized in  table  II  and  figure  11.  Also  included  on  these  displays 
are  the  infinite  media  or  infinite  depth  of  burial  analytical  solu- 
tions described  in  section  I for  the  appropriate  material  param- 
eters, i.e.,  case  A with  NS  at  the  soil-structure  interface. 

An  examination  of  figure  11  indicates  that  for  a depth  of 
burial  of  two  diameters,  the  soil  stresses  are  essentially  the 
same  as  the  infinite  media  soil  stresses  along  the  crown  and 
springline  rays  for  all  of  the  liner  stiffness  considered. 
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Figure  9.  Finite  Element  Mesh  for  SAP4  Static  Analysis,  h/R  = 2.0 
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Table  II 


CYLINDER  NONDIMENS IONAL  THRUST  AND  MOMENT 


FOR 

VARIOUS 

DEPTHS 

OF  BURIAL  - SAP 4 RESULTS 

(R/t  = 15, 

8,6,4) 

\ 

Thrust 

Moment 

h/R 

r/N^ 

0° 

90° 

180° 

0° 

90 

180° 

0.5 

15 

0.514 

1.292 

0.536 

0.01195 

-0.00801 

0.00812 

8 

0.504 

1.314 

0.556 

0.05020 

-0.04420 

0.04430 

6 

0.477 

1.315 

0.549 

0.08430 

-0.08190 

0.08270 

4 

0.428 

1.312 

0.531 

0.13600 

-0.14270 

0.14600 

1.0 

15 

0.540 

1.295 

0.549 

0.00928 

-0.00816 

0.00812 

8 

0.545 

1.331 

0.566 

0.04610 

-0.04450 

0.04450 

6 

0.527 

1.341 

0.556 

0.08290 

-0.08300 

0.08360 

4 

0.490 

1.350 

0.531 

0.14320 

-0.14790 

0.15020 

2.0 

15 

0.554 

1.290 

0.554 

0.00838 

-0.00820 

0.00812 

8 

0.567 

1.340 

0.569 

0.04520 

-0.04500 

0.04490 

6 

0.552 

1.370 

0.556 

0.08460 

-0.08500 

0.08520 

4 

0.519 

1.390 

0.526 

0.15300 

-0.15500 

0.15600 

4.0 

15 

0.555 

1.290 

0.556 

0.00816 

-0.00814 

0.00808 

8 

0.569 

1.350 

0.569 

0.04520 

-0.04510 

0.04500 

6 

0.555 

1.376 

0.555 

0.08600 

-0.08600 

0.08600 

4 

0.521 

1.413 

0.517 

0.15890 

-0.15930 

0.15990 

CO 

15 

0.523 

1.301 

0.523 

0.00867 

-0.00935 

0.00867 

(case 

A)  8 

0.539 

1.364 

0.539 

0.04810 

-0.05050 

0.04810 

6 

0.532 

1.395 

0.532 

0.08970 

-0.09420 

0.08970 

4 

0.512 

1.438 

0.512 

0.16060 

-0.17070 

0.16060 

The  curves  also 

show  that  the 

state 

of  stress  in  the  soil 

immediately 

adjacent  to  the 

structure,  i.e 

. , the 

SSI  stresses,  at  the 

crown  is 

essentially  unaffected  by  the  depth  of  burial  for  the  flexible 
liner  case  (R/ t = 15).  As  the  structural  stiffness  decreases,  the 
SSI  stresses  at  the  crown  are  somewhat  but  not  drastically  reduced. 
There  is  a moderately  pronounced  depth  of  burial  effect  on  the  state 
of  soil  stress  for  points  that  are  slightly  away  from  the  crown. 


The  soil  stresses  along  the  springline  ray  (0  = 90  deg)  are  not 
drastically  affected  by  the  depth  of  burial.  In  fact,  it  can  be 
argued  that  the  springline  ray  soil  stress  differences  between  the 
infinite  media  and  the  SAP4  results  can  be  attributed  to  the  in- 
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25  ksi  (172.4  MPa) 
0.25 


= 3000  ksi  (20,690  MPa) 

= 0,2  HU 


CURVE  h/K 

A * (Case  A-NS) 


Table  II  gives  the  nondimens ional  thrust  and  bending  moment  at 
the  crown,  springline,  and  invert  locations  (9=0,  90,  180  deg) 
for  the  16  SAP4  depth  of  burial  solutions  and  the  infinite  depth 
of  burial  analytical  solution  (case  A) . These  results  show  that 
the  thrusts  and  moments  for  a depth  of  burial  of  two  diameters 
are  generally  well  within  10  percent  of  the  infinite  depth  of 
burial  resulsts.  The  largest  difference,  13  percent,  occurs  for 
the  springline  moment,  and  these  differences  can  again  be  essen- 
tially attributed  to  the  different  solution  schemes  that  were 
employed . 

Using  the  SAP4  results  only,  it  is  observed  that  the  cylinder 
thrusts  are  reduced  as  the  depth  of  burial  is  reduced  for  all  ex- 
cept the  springline  thrust  for  the  most  flexible  cylinder  (R/t  = 15) 
which  increases  very  slightly.  As  the  depth  of  burial  decreases, 
the  crown  moment  increases  for  the  flexible  liner  and  decreases 
for  the  stiffer  liner.  For  the  more  flexible  liners,  the  depth  of 
burial  has  very  little  influence  on  the  springline  and  invert 
moment;  but  for  the  stiffer  liners,  these  moments  slightly  decrease 
(as  much  as  12  percent)  as  the  depth  of  burial  decreases. 

The  overall  depth  of  burial  effect  for  static  loads  on  buried 
structures  appears  to  be  the  negating  influence  or  moderation  of 
the  soil  arching  phenomena.  For  example,  the  soil  arching  around 
a flexible  structure  is  reduced  and  hence  the  loading  is  increased 
as  the  depth  of  burial  decreases.  Conversely,  the  action  of  the 
soil  arching  onto  a stiff  structure  is  also  reduced  as  the  depth 
of  burial  is  reduced. 

For  depths  of  burial  equal  to  one  tunnel  liner  radius  or  more, 
the  primary  SSI  loads  are  essentially  derived  from  modes  0 and  2 
as  evidenced  by  the  Table  II  results.  This  is  apparent  since  the 
crown  and  invert  moment  and  thrust  are  nearly  equal  which  implies 
that  the  uniform  and  ovaling  load  components  dominate.  For  a depth 
of  burial  of  one-half  the  tunnel  radius,  additional  (higher)  modes 
of  the  SSI  loading  are  required  to  accurately  determine  the  load- 
ing; however,  modes  0 and  2 will  still  produce  quite  credible  re- 
sults even  for  this  depth  of  burial  case  for  all  but  the  most 
flexible  liner. 


m 


The  mode  0 and  2 SSI  stresses  can  be  approximated  in  the  fol- 
lowing manner  by  using  the  liner  force  resultants  given  in  table 
II.  First  the  mode  0 and  2 thrust  and  moment  modes  are  calculated 
from  the  crown,  springline,  and  invert  values  in  the  table 
(6  = 0,  90,  180)  as 


^TCR  + TINV^/4  + TSPR/2  1 

I 

To  " ^SPR 
Ho  = ^CR  + MINV^4  + ^SPR^2 
n2  no  mSPR 

Neglecting  the  Mq  solution,  the  desired  SSI  modal  stresses 

(crro  a ~ * rr02^  can  obtained  by  application  of  equations  (lie) 

through  (llg) . 
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SECTION  III 
DYNAMIC  ANALYSES 

* 

The  dynamic  airblast  loading  finite  element  investigations 
and  the  related  processing  of  the  numerical  results  undertaken 
during  this  project  are  described  in  this  section.  'The  fundamen- 
tal purpose  of  these  dynamic  solutions  was  to  see  if  the  important 
aspects  of  the  SSI  loads,  i.e.,  radial  and  shear  stresses,  could 
be  satisfactorily  characterized  by  a few  time  varying  amplitudes 
of  the  modal  components  of  the  SSI  loads.  This  concept  was  pos- 
tulated as  a logical  extension  of  the  technique  discussed  in  the 
static  analyses  in  section  II,  where  it  was  shown  that  modes  0 
and  2 generally  provide  an  accurate  definition  of  the  SSI  loads. 
For  the  static  analyses,  both  the  tunnel  liner  and  the  soil  media 
were  idealized  as  elastic  materials  and  the  results  are  propor- 
tional  to  the  loading.  Since  the  dynamic  analyses  employ  a non- 
linear soil  material,  the  solutions  are  sensitive  to  the  magnitude 
of  the  loading  and  therefore  both  high  and  low  values  of  peak  air- 
blast  overpressure  have  been  included  in  these  analyses.  These 
dynamic  investigations  use  the  same  range  of  tunnel  liner  radius 
to  thickness  ratios  and  depths  of  burial  that  were  used  in  the 
static  analyses. 

A plane  strain  solution  procedure  has  been  used;  in  general, 
the  soil  medium  has  been  modeled  as  a nonlinear  material  and  the 
buried  cylinder  liner  was  representative  of  a linear  elastic  con- 
crete material  with  a 13.12  ft  (4.0  m)  inside  diameter  (see  fig- 
ure 12) . The  dynamic  loading  is  assumed  to  originate  from  a 
nuclear  weapon  airblast  overpressure.  The  passage  of  the  airblast 
shock  front  over  the  ground  surface  has  not  been  included  in 
these  in  'tigations,  and  the  airblast  loading  is  applied  uni- 
formly to  the  entire  ground  surface.  Referring  to  figure  12  the 
airblast  loading  is  represented  by  a time  varying  pressure  acting 
uniformly  over  the  ground  surface.  Table  III  presents  the  time 
duration  of  an  equivalent  triangular  pulse  of  intense  airblast 
overpressures  for  a reasonable  range  of  weapon  yields,  0.125, 
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Figure  12.  Dynamic  Analysis  Parameters. 
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1.0,  and  8.0  MT  (0.53xl09,  4.2xl09,  and  33.6xl09  HJ)  and  peak 
overpressures,  pQ,  of  500,  1000,  and  1500  psi  (3.448,  6.395,  and 
10.343  MPa) . 


Table  III 

APPROXIMATE  EQUIVALENT  TRIANGULAR  PULSE  AIRBLAST 
OVERPRESSURE  DURATION  (msec) 


Overpressure,  p 

Weapon  Yield, 

MT  (MJ) 

psi  (MPa) 

0.125  (0.525xl09)  1.0 

(4.2xl09) 

8.0  (33.6xl09) 

500  (3.448) 

too 

C50 

10 

20 

40 

25 

50 

100 

1000  (6.895) 

t"00 

50 

8 

15 

30 

15 

30 

60 

1500  (10.343) 

coo 

C50 

5 

10 

20 

10 

20 

40 

The  duration  parameters  were  derived  from  reference  49  for  two 
triangular  idealizations: 

t^p,  derived  by  using  the  slope  of  the  triangular  pulse 
as  equal  to  the  initial  decay  slope  of  the  actual 
pressure  time  curve. 

t^g,  derived  by  assuming  that  the  triangular  pulse  passes 
through  the  actual  curve  at  0.5  pQ . 

The  primary  objective  of  these  dynamic  solutions  was  to  obtain 
SSI  loads  acting  on  the  structure  and  internal  liner  forces  for 
a reasonable  variation  of  the  significant  soil  and  structure  param- 
eters . 

1.  RUN  MATRIX  OF  DYNAMIC  CASES 

A desirable  matrix  of  solution  parameters  that  might  be  se- 
lected consistent  with  the  static  investigations  could  conceiv- 
ably include  three  or  four  variations  for  each  fundamental  param- 
eter. For  example,  consider  a run  matrix  that  includes  these  five 
parameters : 
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© Structural  Stiffness  (elastic  concrete) 

R/t  = 15,  8,  6,  4 

e Depth  of  Burial 

h/R  = 0.5,  1.0,  2.0,  4.0 

® Peak  Airblast  Overpressure 

Po  = 500,  1000,  1500  psi  (3.45,  6.90,  10.34  MPa) 

e Airblast  Loading  Duration 

td  = 0.015,  0.030,  0.060  sec 

• Soil  Media  (see  figure  13) 

Sand,  Cemented  Sand,  Clay. 

Even  this  rather  restrictive  set  of  parameters  would  result  in 
4x4x3x3x3  = 432  cases  if  all  combinations  were  taken.  Obviously 
this  number  of  cases  would  have  resulted  in  an  effort  that  was 
well  beyond  the  scope  of  work  for  this  research  project.  To  study 
the  effects  of  the  above  parameter  ranges  on  the  dynamic  response 
of  the  tunnel  liner,  a series  of  cases  were  made  up  using  the  ex- 
treme values  of  the  above  nonsoil  parameters: 


R/t  = 15  and  4 
h/4  = 0.5  and  4.0 

Po  = 500  and  1500  psi  (3.45  and  10.34  MPa) 
t^  = 0.015  and  0.060  sec. 

Parameters  associated  with  two  soil  media  (sand  and  cemented  sand) 
were  also  used  to  define  the  dynamic  run  matrix.  All  combinations 
of  these  parameters  result  in  2x2x2x2x2  = 32  cases.  Four  addi- 
tional cases  were  also  defined  which  involve  elastic  soil  media 
and  nonlinear  structure  parameters.  The  resulting  matrix  of  solu- 
tion cases  defining  each  of  the  36  dynamic  finite  element  solutions 
is  given  in  table  IV. 
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Table  IV 

RUN  MATRIX  FOR  (SAMSON)  DYNAMIC  SOLUTIONS 


Case 

R/t 

h/R 

Po 

psi 

Cd 

msec 

Soil3 

Case 

R/t 

h/R 

Po 

psi 

Cd 

msec 

Soil3 

1 

4 

0.5 

500 

15 

A 

21 

4 

4 

500 

15 

B 

2 

4 

0.5 

1500 

22 

4 

4 

1500 

3 

4 

0.5 

500 

60 

23 

4 

4 

500 

60 

4 

4 

0.5 

1500 

24 

4 

4 

1500 

5 

4 

0.5 

500 

15 

B 

25 

15 

4 

500 

15 

A 

6 

4 

0.5 

1500 

26 

15 

4 

1500 

7 

4 

0.5 

500 

60 

27 

15 

4 

500 

60 

8 

4 

0.5 

1500 

28 

15 

4 

1500 

9 

15 

0.5 

500 

15 

A 

29 

15 

4 

500 

15 

B 

10 

15 

0.5 

1500 

30 

15 

4 

1500 

11 

15 

0.5 

500 

60 

31 

15 

4 

500 

60 

12 

15 

0.5 

1500 

32 

15 

4 

1500 

13 

15 

0.5 

500 

15 

B 

33c 

15 

0.5 

1500 

15 

Elastic1 

14 

15 

0.5 

1500 

34C 

15 

0.5 

1500 

15 

B 1 

15 

15 

0.5 

500 

60 

35 

4 

0.5 

500 

15 

Elastic 

16 

15 

0.5 

1500 

36c 

4 

0.5 

500 

15 

B 

17 

4 

4 

500 

15 

A 

18 

4 

4 

1500 

19 

4 

4 

500 

60 

20 

4 

4 

1500 

aSoil  A - Sand;  Soil  B - Cemented  Sand 
^Elastic  Soil;  E = 25  ksi,  V = 0.25 


CNonlinear  Liner, 

1 

percent 

Steel  R/C 

Note: 

Cases  33, 

34 

same  as 

14  except  as 

indicated. 

Cases  35, 

36 

same  as 

5 except  as 

indicated . 

1 psi 

= 6.895  kPa 
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2.  DYNAMIC  FINITE  ELEMENT  SOLUTION  ANALYSIS 


A modified  version  of  the  SAMSON  computer  code  (Ref.  48)  was 
used  to  obtain  solutions  to  the  36  cases  defined  in  table  IV. 

This  computer  code  is  a finite  element  code  for  the  two-dimesional 
dynamic  stress  analysis  of  plane  and  axisymmetric  media.  It  is 
based  on  small  strain  theory  and  is  particularly  suited  for  han - 
dling  problems  with  nonlinear  material  properties  and  a large  num- 
ber of  degrees  of  freedom.  Although  it  was  designed  especially 
for  the  investigation  of  SSI  problems,  it  can  be  used  for  more  gen- 
eral applications. 

For  nonlinear  problems  a direct  method  of  solution  is  employed 
without  explicit  formulation  of  a stiffness  matrix.  This  method 
has  several  advantages  for  large-scale  nonlinear  transient  prob- 
lems. By  avoiding  construction  of  the  stiffness  matrix,  consider- 
able computer  storage  is  saved,  and  the  solution  of  much  larger 
problems  is  permitted.  Moreover,  it  is  also  computationally  more 
efficient  than  other  commonly  used  solution  methods  in  nonlinear 
analysis.  For  linear  elastic  problems,  the  normal  stiffness  matrix 
method  of  solution  can  be  optionally  used.  The  following  fea- 
tures of  the  SAMSON  code  were  used  to  obtain  the  solutions  to  the 
36  dynamic  cases . 

a.  Automatic  Generation  of  Mesh  Data 

The  SAMSON  code  has  a mesh  generation  capability  which  allows 
the  user  to  skip  node  and  element  cards  in  the  data  deck.  The 
code  fills  in  missing  node  cards  by  assuming  uniform  spacing  be- 
tween the  nodes  for  which  data  are  provided.  Missing  element 
cards  are  generated  by  assuming  that  the  missing  node  numbers  may 
be  derived  from  the  node  numbers  of  the  previous  element. 

b.  Pressure  and  Force  Loadings 

Arbitrary  time  variations  in  pressure  or  forces  can  be  speci- 
fied at  nodes.  The  forces  represented  by  pressures  at  nodes  are 
computed  assuming  a linear  variation  in  pressure  between  nodes 
along  a user-defined  surface. 
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Materials 


c . 

Two  types  of  material  constitutive  models  were  used  in  the 
analysis,  elastic  and  nonlinear.  The  elastic  model  incorporates 
a linear  elastic  material  which  was  used  to  model  the  concrete 
cylindrical  material  for  the  majority  of  the  cases. 

The  important  parameters  that  define  the  elastic  cylinder  are: 

Ec  = 3000  ksi  (20,700  MPa) 

v = 0.2 

C o 

Yc  = 150  pcf  (2.403  Mg/mJ) , density 

The  nonlinear  soil  medium  was  modeled  by  using  the  soil/rock 
material  model  (R0CKMT3)  described  in  reference  55.  Input  to 
this  material  model  includes  the  density,  loading  and  unloading 
bulk  moduli  and  Poisson's  ratios,  and  a failure  surface  description. 
The  bulk  moduli  are  specified  by  using  a piecewise  linear  variation 
to  approximate  the  mean  normal  stress  (a)  and  the  volumetric  strain 
(e)  as  shown  in  figure  13a.  Tension  stress  cutoffs  can  also  be  in- 
cluded in  the  material  response.  The  yield  or  failure  surface  is 
defined  by  a piecewise  linear  approximation  of  the  /L)  versus  a 
curve  as  indicated  in  figure  13b.  Note  that  only  a very  simple 
"single"  linear  variation  was  used  to  approximate  the  yield  sur- 
face for  the  soil  models  of  this  investigation.  The  variables 
used  in  the  R0CKMT3  material  model  are  defined  as: 

a = (Oi  + a2  + Og)/3  (12a) 

e = ei  + e2  + e3  (12b) 

J2  = \ Si  + S2  + S3  I (12c) 

where  a ^ and  are  principal  components  of  stress  and  strain 

S.  = a.  - a,  deviatoric  stresses  (12d) 

ii’ 

The  $lope  of  each  piecewise  linear  material  response  characteristic 
curve  (Aa/Ae)  is  the  bulk  modulus  (K)  of  the  material  for  each 
state  of  loading  or  unloading.  Figure  14  gives  dynamic  uniaxial 
test  results  (Ref.  56)  that  have  been  used  to  define  simplified 
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R0CKMT3  stress-strain  characteristics  for  both  of  the  sand  mate- 
rial models  (A  and  B)  used  in  the  dynamic  finite  element  analyses 
It  has  been  tacitly  assumed  that  the  unloading/reloading  path  is 
parallel  to  the  highest  or  last  slope  of  the  curve,  i.e.,  for 
a^<-1.5  ksi  (-10.34  MPa).  The  confined  compression  moduli  (M) 
for  each  of  the  three  linear  rates  used  in  the  approximation  are 
given  on  the  uniaxial  data  curve.  Poisson's  ratio  for  both  sand 
materials  (A  and  B)  have  been  taken  as  v = 0.25.  Thus  an  elastic 
biaxiality  ratio  of  lateral  to  axial  stress  is  k = v/ (1-v)  = 1/3. 
The  uniaxial  approximation  data  of  figure  14  were  used  to  obtain 
the  figure  13a  curves.  For  example,  the  (a,  e)  point  at  the  end 
of  the  first  linear  portion  is  obtained  by  application  of  equa- 
tion (12a)  and  (12b)  where 


There  results 


o-L  = -0.5 

°2  = °3  = ka1  = -0.5/3 
e1  = -0.0167 
e2  = c3  = 0 

a = (l+2k)0l/3 
= -0.279 


e 


el 

-0.0167 


(13a) 


(13b) 


Thus  the  volumetric  strains  at  the  break  points  on  the  piecewise 
linear  o versus  e curve  are  the  same  as  the  axial  strains  (e^)  for 
the  uniaxial  approximation.  The  bulk  modulus  (K)  is  computed  from 
the  modulus  of  confined  compression  (M)  from  equations  (13a, 13b). 
There  results 


K = M(l+2k)/3 


5 

1 


M 


(13c) 
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Referring  to  figure  15,  the  relationship  between  the 
versus  o failure  surface  and  the  angle  of  internal  friction  (#) 
and  cohesion  (c)  of  the  material  is  now  determined.  Consider 
a triaxial  shear  test  state  of  stress  (figure  15a)  where  the  axial 
pressure  is  p-^  and  the  hydrostatic  lateral  pressures  are  P2  = P3 
where  p^  = -o^.  The  stress  state  on  any  plane  from  the  vertical 
plane  (9)  is 


2 2 

p = p^  cos  9 + P3  sin  0 


(14a) 


t = (p-^-p^)  cosQ  sinQ 
= 0.5  (P3-P3)  sin29 


(14b) 


The  Coulomb-Mohr  law  for  shear  failure  (figure  15c)  can  be  written 


t = c + p tan^ 


(14c) 


where  c and  0 are  determined  from  the  envelope  of  the  Mohr  diagram 
circles  drawn  from  the  triaxial  shear  tests.  Substituting  equa- 
tions (14a)  and  (14b)  into  (14c)  and  rearranging  terms  yields 


Pi  " 


c+p3(sin0  cos9  + sin  9 tan£) 


sin9  cos9  - cos  9 tan^ 


(15a) 


The  critical  plane  (0cr)  is  found  when  p-^  is  a minimum  to  cause 
failure;  and  p-^  is  a minimum  when  the  denominator  of  equation  (15a) 
is  a maximum.  This  occurs  when  the  rate  of  change  of  the.  denomina- 
tor with  respect  to  9 is  zero.  Performing  the  indicated  operation 
on  the  denominator  results  in  the  well  known  solution  of 


9cr  * f+«° 


(15b) 


From  elementary  trigonometry 


tan$  = cot(90-$) 


= -cot (90+0) 


= -cot  20, 


(15c) 
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(a)  TR I AXIAL  SHEAR  TEST 


(b)  STRESS  TRANSFORMATION 


Using  equation  (15b)  for  0 in  equation  (15a)  and  also  employing 

equation  (15c)  and  multiplying  the  numerator  and  denominator  of 

the  equation  (15a)  by  2tan0  , it  can  be  shown,  by  using  the 

cr  2 

trigonometric  identity  cot  20  = (1-tan  0)/(2  tan0) , that 

p^  = 2c  tan(2  + 45)  + p3  tan2(^+45)  (15d) 

For  triaxial  tests,  the  principal  stresses  are  the  axial  (o-^  = -p^) 
and  the  hydrostatic  lateral  stress  (02  = 03  = -P2  = -p3) ; there- 
fore , 

o = - 3 (pi  + 2p3)  (16a) 


S 


1 


- a 

(°r°  3) 


So  = s,  - 


o3  -o 


t (0ra3) 


= - — s 
2 1 


J2 


1 

" 7 
. 1 


2 2 
Sj  + 2S^ 


(0r°3) 


SJ'2  - (P1"P3)/’/^  (16b) 


rr 


Solving  for  from  equation  (16a)  and  then  using  this  result  and 
equation  (16d)  in  (16b)  gives 


o 4- 


(17a) 


l 


Referring  to  figure  15d  it  is  observed  that  equation  (17a)  is  of  the 
form  y = mx  + b where  b is  the  intercept  of  the  straight  line  with 
the  /j^,  axis  and  m is  the  slope  of  the  line.  From  equation  (17a) 


m 


= /5 


1-t 

2+t‘ 


2 \ 


(17b) 

j 


b = 2/3  — ^ c 

2+tz 


(17c) 


For  both  of  the  sand  materials  considered  (figure  13b) , the  slope 
of  the  failure  surface  is  m = -0.85.  Using  this  value  in  equation 
(17b)  and  solving  for  t and  then  using  equation  (16c)  to  calculate 
the  conventional  angle  of  internal  friction  ($)  yields 

<t>  - 36.234°,  Soil  materials  A and  B 


Since  b = 0 for  the  sand  material  (A) , the  cohesion  (c)  is  zero  by 
equation  (17c).  However  b = 0.1  ksi  (0.6895  MPa)  for  the  cemented 
sand  material  (B)  and  this  results  in  a cohesion  magnitude  of 

c = 0.0862  ksi  (0.5943  MPa),  soil  material  B 

c = 0,  soil  material  A 

3 

A soil  density  of  y = 109.6  pcf  (1.756  Mg/m  ) was  used  for  both 
materials . 

The  nonlinear  reinforced  concrete  cylinder  used  for 
cases  34  and  36  was  modeled  by  also  using  the  R0CKMT3  material 
model.  A fiber  reinforced  concrete  material  with  1 percent  steel 
(Ref.  57)  has  been  selected  for  this  liner  idealization.  Figure  16 
gives  the  material  characteristics  for  this  model. 

d.  Damping 

The  SAMSON  Code  provides  an  option  for  incorporating  both  linear 
damping,  proportional  to  the  strain  rate  in  an  element,  and  nonlin- 
ear damping,  proportional  to  the  square  of  the  strain  rate.  Both 
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Figure  16.  Nonlinear  R/C  Liner  (1%  Steel)  R0CKMT3  Parameters 


types  of  damping  consist  of  a hydrostatic  component  and  a devia- 
toric  component  which  may  be  applied  separately  or  in  combination. 
The  following  damping  parameters  have  been  found  to  yield  good 
results  for  buried  structure  dynamic  analyses  by  AFWL  personnel. 

(1)  Linear  Coefficients 

b^  = 0.04,  hydrostatic 
b2  = 0.02,  deviatoric 

(2)  Nonlinear  Coefficients 

b^  = b^  = 2.0,  hydrostatic  and  deviatoric 

Coefficients  b^  and  b2  are  approximately  the  fractions  of  criti- 
cal damping. 

e.  Finite  Element  Model 

The  finite  element  models  of  the  soil  media  for  the  two  depths 
of  burial  considered  (h/R  = 0.5  and  4.0)  are  shown  in  figures  17 
and  18.  These  meshes  are  very  similar  to  the  static  analysis  meshes 
(figures  7 and  10)  except  more  refinement  of  the  gridwork  has  been 
used  in  certain  portions  and  the  entire  mesh  is  composed  of  tri- 
angular and  rectangular  elements.  The  latter  is  required  since  the 
SAMSON  code  does  not  have  a general  quadrilateral  finite  element 
capability.  The  added  mesh  refinement  was  incorporated,  especially 
between  the  ground  surface  and  the  tunnel  liner,  so  that  the  air- 


blast  induced  stress  wave  traversing  through  the  soil  would  not  be 
unduly  modified  by  abrupt  element  size  discontinuities. 

The  left  and  right  boundaries  of  each  mesh  are  free  to  move 
vertically,  and  the  lower  boundary  is  free  to  move  horizontally. 
The  lower  boundary  is  suf f iciently  removed  from  the  cylindrical 
liner  so  that  reflected  waves  from  this  boundary  do  not  reach  the 
structure  in  the  time  of  interest  for  the  short  duration  airblast 
loads  (0.015  sec).  For  the  long  duration  loads  (0.060  sec)  addi- 
tional elements  were  added  to  the  lower  portion  of  the  mesh  to 
extend  the  lower  boundary  by  two  outside  tunnel  liner  diameters  or 
by  4.0  R. 
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Mesh  4DS  (liner  not  shown) 

354  Nodes 

608  Elements 

Triangular  finite  elements  have  also  been  used  to  model  each 
‘ of  the  tunnel  liners.  The  outer  nodes  of  the  liner  mesh  were  at- 

tached to  the  soil  mesh  so  that  no  relative  motion  was  permitted 
at  the  soil-structure  interface.  In  an  effort  to  comprehensively 
model  both  the  membrane  and  the  bending  deformations  of  the  tunnel, 
the  liner  thickness  w r subdivided  by  four  layers  of  equal  thick- 
ness The  circumference  was  subdivided  by  rays  at  7.5  deg  incre- 
ments, which  yields  24  equal  divisions  for  the  180  deg  half-model 
of  the  liner.  This  results  in  4x24  = 96  cells  requiring  two  tri- 
angular elements  per  cell  for  a total  of  2x96  = 192  elements  to 
model  each  tunnel  liner.  Figure  19  shows  the  resulting  finite 
element  mesh  for  the  most  flexible  liner  (R/ t = 15)  and  a portion 
of  the  soil  media  mesh.  Each  of  the  24  circumferential  liner  sub- 
divisions has  been  numbered  as  indicated  on  this  figure  for  future 
reference  in  the  numerical  solutions . The  finite  element  mesh  for 
the  stiff  tunnel  liner  (R/t  = 4)  is  defined  in  the  same  way  as  the 
flexible  liner  (i.e.,  four  subdivisions  through  the  thickness  and 
24  circumferential  segments  requiring  192  triangular  elements)  with 
the  obvious  difference  that  the  node  points  are  further  apart, 
through  the  thickness,  to  model  the  increased  liner  thickness. 

As  previously  noted,  an  inside  radius  of  the  tunnel  liner  of 
6.56  ft  (2  m)  was  maintained  for  all  of  the  SAMSON  numerical  solu- 
tions. This  was  implemented  by  providing  basic  mesh  node  point 
data  so  that  the  radius  to  the  soil-structure  interface  was  1.0  and 
then  applying  an  appropriate  scale  factor  to  the  coordinates,  with- 
in the  SAMSON  code,  so  that  the  inside  radius  of  the  liner  was 
6.56  ft  (2  m) . 

f.  SAMSON  Modifications 

The  output  element  stresses  in  the  computer  code  are  related  to 
the  vertical  (z)  horizontal  (x)  coordinate  system.  The  SSI  stresses 
and  the  tunnel  liner  internal  force  resultants  are  the  desired  solu- 
tion results.  To  facilitate  this,  program  modifications  were  made 
such  that  the  element  stresses  are  referred  to  the  radial  (r)  and 
circumferential  (9)  coordinates  for  the  elements  of  interest,  i.e., 
the  192  tunnel  liner  elements  and  the  48  triangular  soil  elements  in 

i 
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Figure  19.  Tunnel  Liner  Mesh  for  SAMSON 
Dynamic  Analysis  (R/t  = 15) 
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the  24  cells  immediately  adjacent  to  the  structure.  The  SSI 
stresses  (or,trQ)and  also  the  circumferential  stress  (°q)  were  ob- 
tained at  each  of  the  24  reference  points  (figure  19)  by  computing 
the  average  stress  in  the  two  triangular  soil  elements  at  each 
reference  point. 

The  tunnel  liner  force  resultants  (M,T,V)  at  the  24  reference 
points  were  computed  by  first  converting  the  liner  element  stresses 
to  the  r,9  system  and  then  obtaining  the  stress  in  each  of  the 
four  cells,  through  the  liner  thickness,  as  the  average  stress  of 
the  two  triangular  elements  in  the  cell.  Referring  to  figure  19 
denote  the  hoop  stress,  Oq,  for  each  of  the  four  cells  at  a given 
reference  point,  I,  as  ai*a2,03>a4  where  the  subscripts  1 through  4 
correspond  to  the  inner  through  outer  cells,  respectively.  Then 
the  thrust  and  moment  at  the  reference  point  can  be  computed  as 


T = - (oj.  + a 2 + a3  + o^)  t/4  (18a) 

M = (3^  +a2  - a3  - 3<>4)  t2/32  (18b) 

In  a similar  manner,  the  liner  shear  is 


V = - (t^ + i2  + T3  + T4)  t/4 


(18c) 


where 


is  the  average  shear  stress  (x^g)  in  the  ith  cell. 


Nondimens ional  liner  force  resultants  and  also  nondimens ional 
SSI  stresses  were  computed  at  each  of  the  24  reference  points  at 
uniform  time  increments  (4tQut)  during  the  SAMSON  solution.  These 
results  were  written  on  a magnetic  tape  for  each  of  the  36  cases 
for  use  as  input  to  the  postprocessor  codes.  For  the  more  flexi- 
ble liner  (R/t  = 15)  cases  the  results  were  stored  on  tape  at 


At  . = 0.1  msec  and  for  the  liner  with  R/t  = 4, 
out 

3.  RESULTS  OF  DYNAMIC  (SAMSON)  SOLUTIONS 


At  _ = 0.075  msec, 
out 


A postprocessing  code  was  written  to  facilitate  the  processing 
and  reporting  of  the  SAMSON  numerical  solution  results.  Both  tab- 
ular printout  and  plotted  time  history  curves  were  obtained  for 
all  36  cases.  All  of  the  results  are  nondimensional,  and  any 


future  reference  to  the  liner  internal  forces  and  SSI  stresses 
refers  to  the  nondimensional  values  unless  otherwise  stated. 

Nondimensional  Liner  Force  Resultants 

Thrust , T/p  R — > T 
o o 

Moment,  M/pQR  — > M 

Shear,  V/pQR  — v V 

I ■■ 

Nondimensional  SSI  Stresses 
Radial,  ar/pQ  — *■  STR-R 
Circumferential,  ag/p0— >STR-T 
Shear,  TrQ/po TAU-RT 

The  majority  of  the  tabulated  and  plotted  results  for  the  36 
cases  defined  in  table  IV  are  presented  in  volume  u 0f  this  report. 

The  results  are  grouped  for  each  case,  and  the  tables  are  presented 
first  followed  by  the  GOULD  plotted  curves.  Each  table  and  graph 
in  volume  II  contains  information  to  identify  the  data  and  the  case 
number. 

The  results  for  case  5 are  presented  at  this  point  to  illus- 
trate the  kind  of  information  that  has  been  obtained  for  all  cases. 

a.  Peak  Response  Results 

The  minimum  and  maximum  values  of  the  nondimensional  liner  in- 
ternal force  resultants  and  the  SSI  nondimensional  stresses  that 
occurred  at  each  of  the  24  reference  points  during  the  complete 
solution  history  are  presented  in  table  V.  Reference  points  1,  12, 
and  24  are  denoted  as  the  crown,  springline  and  invert  points; 
however,  they  actually  refer  to  the  center  of  the  segment  that  is 
closest  to  the  9 = 0,  90,  and  180  deg  points. 

Figures  20a  through  20f  give  curves  of  the  maximum  and  minimum 
results  of  table  V.  In  these  figures,  the  abscissa  title  Liner-Pt. 
is  the  same  as  the  24  reference  points  (see  figure  19)  around  the 
tunnel  liner. 

b.  Modal  History  Results 

As  previously  discussed,  the  structural  system  and  the  loading 
environment  for  the  cases  considered  in  this  project  are  symmet- 
ric about  a vertical  plane  through  the  longitudinal  axis  of  the 
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CASE  5 

Figure  20a.  Minimum  and  Maximum  Nondimens ional  Thrust 
Variation  around  Liner,  Case  5 
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CASE  5 

Figure  20c.  Minimum  and  Maximum  Nondimens ional  Shear 
Variation  around  Liner,  Case  5 
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Figure  20e.  Minimum  and  Maximum  SSI  Nondimensional 
Circumferential  Stress  Variation  around  Liner,  Case  5 
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Figure  20f.  Minimum  and  Maximum  SSI  Nondimens ional 
Shear  Stress  Variation  around  Liner,  Case  5 


tunnel.  Therefore,  the  radial  and  circumferential  motions  of  the 
structure  are  symmetric  and  antisymmetric,  respectively,  with  re- 
spect to  the  plane  of  symmetry.  The  liner  thrust  and  moment  and 
the  SSI  radial  and  circumferential  stresses  are  symmetric,  and 
the  liner  shear  and  SSI  shear  stresses  are  antisymmetric.  There- 
fore, each  response  variable  of  interest  can  be  represented  by  an 
appropriate  Fourier  sine  or  cosine  series.  Referring  to  figure  21 
these  responses  can  be  represented  as  follows. 


(1)  Symmetric  Response  (S) 


S(t,9)  = 


s„(t) 


y Sn(t)  cos  n9 
n=l 


(19a) 


where  Sn  e Sn(t)  are  the  time  dependent  modal  amplitude  histories. 
Variables  that  can  be  represented  in  this  manner  are  the  liner 
thrust  and  moment  and  the  soil  media  radial  and  circumferential 
stresses  at  the  SSI  interface.  The  modal  amplitudes  (Sn)  are  eval- 
uated by  multiplying  equation  (19a)  by  cosn9  and  integrating  from 
9 = 0 to  tt,  and  results  in 


S(t,9)  cos  n9  d9 


tt 


cos  n9  d9  n>9 


Therefore 


S =2f 
n it J 


S(t, 9)  cos  n9  d9 


(19b) 


At  any  time  step  t in  the  solution,  the  variable  S(t,9)  is 
known  at  the  24  discrete  reference  points  around  the  tunnel  liner. 
Denoting  these  magnitudes  at  the  24  reference  points  (9^)  as 
S(t,9^),  the  right-hand  side  of  equation  (19b)  can  be  numerically 
integrated,  and  the  modal  amplitude  is 


S = — A9  Y'  S(t,9.)  cos  n9.  , 
n tt  £_  i i ’ 

i=l 


n > 0 


(19c) 
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Figure  21 


iiuLdtion  tor  nuua  i ituaes 

from  Dynamic  (SAMSON)  Results 
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The  zero  mode  amplitude  is  given  by 


V 


24 

S = - A9  Y S (t , 9 . ) (19d) 

O TT  L-.  1 

i=l 

(2)  Antisymmetric  Response  (A) 

00 

A(t , 9)  = y'  An  (t)  sin  n6  (20a) 

n=l 

where  An  = An(t)  are  the  time  dependent  modal  amplitude  histories. 
Variables  that  can  be  represented  in  this  manner  are  the  liner 
shear  force  and  the  soil  media  shear  stresses  at  the  SSI  inter- 
faces. The  modal  amplitudes  (An)  are  evaluated  by  multiplying 
equation  (20a)  by  sinnG  and  integrating  from  9 = 0 to  n , and 
results  in 

TT 

A n = vf  sin  n0  (20b) 

o 

Following  the  numerical  integration  solution  process  used  in  the 
symmetric  response  analysis,  the  modal  amplitude  is  obtained  from 

24 

A = — A9  y A(t,9.)  sin  n9.  , n>0  (20c) 

n it  4—  l i 

i=l 

The  above  techniques,  equations  (19c)  and  (20c),  were  used 
to  convert  the  dynamic  results  obtained  from  the  SAMSON  code  to 
modal  amplitude  histories. 

(3)  Symmetric  Modal  Histories 

Modal  histories  were  computed  for  modes  0 through  6 for  the 
following  symmetric  response  variables.  The  modal  history  curves 
for  each  mode  are  given  in  the  indicated  figures. 

Nondimensional  Liner  Thrust,  figure  22  (Sn^-Tn/pQR) 

o 

Nondimensional  Liner  Moment,  figure  23  (Sj^-M^p^  ) 

Nondimensional  SSI  Radial  Stress,  figure  24  (Sn^-arn/pQ) 


73 


THRUST r MODE  1 
-1.00  .00  1.00 


(4)  Antisymmetric  Modal  Histories 

Since  the  zero  mode  does  not  exist  for  the  antisymmetric  modes, 
the  modal  histories  were  computed  for  modes  1 through  6 for  the 
following  variables. 

Nondimens ional  Liner  Shear,  figure  25  (An<_Vn/pQR) 

Nondimensional  SSI  Shear  Stress,  figure  26  (An^—  TrQn/P0) 

Modal  amplitude  histories  for  the  circumferential  soil  stresses 
were  not  determined  since  these  stresses  do  not  act  as  loads  on 
the  liner  structure  and  are  only  of  academic  interest. 

Table  VI  gives  the  minimum  and  maximum  values  of  the  modal 
amplitude  that  occur  during  the  time  history  for  each  mode.  Since 
these  peaks  do  not  occur  at  the  same  time  for  all  modes,  the  maxi- 
mum and  minimum  values  of  response  at  a particular  point  (9^)  can- 
not be  obtained  by  simply  summing  the  angular  modification  of  the 

peak  values,  i.e.,  £S  cos  n0.  or  £ A sin  n9..  Even  if  the  peaks 
do  occur  at  the  same  time,  the  maximum  value  of  9^  is  in  general 
not  obtained  from  the  sum  of  the  maximum  modal  amplitudes  since 
the  sin  n0^  or  cos  n9^  will  change  the  effect  from  maximum  at  the 
crown  (9  = 0)  to  minimum  (or  negative)  at  9^  depending  on  the 
sign  of  the  appropriate  trigonometric  function.  For  example,  if 
all  peak  amplitudes  (minimum  and  maximum)  occurred  at  the  same 
time  for  a symmetric  variable,  the  maximum  value  of  the  variable 
at  say  9 = 90  deg  would  be  obtained  (for  the  modes  0 through  6)  as 

6 5 

Max.  S(9  = 90)  = ^ Max.  Sn  + ^ Min.  Sn 

n=  0 , 2 , n= 1 , 3 , 

c.  Crown,  Springline,  Invert  Results 

Time  history  plots  of  the  liner  nondimensional  moment  and 
thrust  and  the  nondimensional  SSI  radial  stress  computed  during 
the  dynamic  SAMSON  analysis  are  presented  for  three  reference 
points  on  the  liner.  The  reference  points  and  figures  containing 
the  history  plots  are: 

Crown  - Ref.  Pt.  1,  9^  = 3.75  deg,  figure  27 
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Table  VI 

MINIMUM  AND  MAXIMUM  FOURIER  MODAL  AMPLITUDES  OF  NONDIMENS IONAL 
LINER  FORCE  RESULTANTS  (T.M.V)  AND  SSI  STRESSES  (a  . t ) CASE  5 
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CASE  5 

Figure  27a.  Nondimensional  Crovm  Thrust  History,  Case  5 


Springline  - Ref.  Pt.  12,  9^  = 86.25  deg,  figure  28 

Invert  - Ref.  Pt.  24,  4 = 176.25  deg,  figure  29 

Hereafter,  any  allusion  to  the  crown,  springline,  and  invert  points 
will  actually  correspond  to  the  above  identified  reference  points 
unless  otherwise  indicated. 

The  final  bit  of  processing  that  was  performed  on  the  SAMSON 
results  was  directed  at  indicating  which  modes  contribute  the  most 
(or  are  the  most  important)  when  trying  to  approximate  the  ele- 
mentary SAMSON  history  results  by  an  appropriate  combination  of  the 
modal  histories.  Table  VII  gives  the  results  of  this  investiga- 
tion, and  a discussion  of  the  parameters  given  therein  follows. 
Processing  was  performed  on  the  nondimens ional  liner  moment  and 
thrust,  and  the  SSI  radial  stress  at  three  reference  points,  i.e., 
crown,  springline  and  invert. 

The  results  presented  in  table  VII  are  described  by  discus- 
sing the  Crown  Thrust  results.  First  consider  the  last  row  of  the 
crown  thrust  results  which  is  labeled  MAX.  The  last  column  in 
this  row  is  labeled  INPUT  VARIABLE  with  a magnitude  of  0.60422. 

This  is  the  maximum  value  of  the  reference  point  1 nondimens ional 
liner  thrust  that  was  obtained  during  the  SAMSON  solution  and  is 
the  same  value  given  in  table  V,  also  see  figure  27a.  If  only  the 
mode  0 nondimens ional  thrust  history  is  used  (see  figure  22a)  that 
was  obtained  by  the  process  described  in  the  derivation  of  equation 
(19d) , the  result  is  a single  mode  (mode  0)  approximation  to  the 
crown  thrust.  The  maximum  nondimensional  thrust  for  this  single 
mode  approximation,  which  for  this  case  is  simply  the  mode  0 max 
given  in  table  VI,  is  0.54571.  The  ratio  of  the  modal  history 
approximation  maximum  to  the  SAMSON  history  maximum  is  0.54571/ 
0.60422  = 0.903163  which  is  the  value  given  in  the  Mode  0 column 
of  table  VII.  The  remaining  columns  for  the  MAX  row  are  similarly 
determined  except  the  approximation  to  the  thrust  is  determined 
by  including  the  modes  0-n  as  indicated  by  the  column  headings. 
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Table  VII 

COMPARISON  OF  MODAL  HISTORY  APPROXIMATION  WITH  SAMSON  SOLUTION  HISTORY  (INPUT  VARIABLE) 
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For  example,  if  modes  0 through  5 are  used  to  determine  the  approx- 
imation to  the  crown  thrust  (denoted  by  Tq  5>  the  maximum  crown 
thrust  derived  from  these  modes  would  be  obtained  from 

Max  Tq  5 = Max  J Tn^ti^cosn0if ’ i=1>L 

ln=0  J 

where  T (t.)  is  the  discrete  value  of  the  mode  n nondimensional 
n 1 

thrust  at  the  ith  time  step  and  there  are  a total  of  L time  steps. 
This  maximum  divided  by  the  SAMSON  max  (0.60422)  would  give  the 
value  given  in  the  Mode  0-5  column,  which  is  0.98075.  If  the  SAMSON 
peak  value,  given  in  the  INPUT  VARIABLE  column,  is  zero  (0)  then  the 
peak  value  obtained  from  the  modal  history  approximation  is  not  nor- 
malized since  this  would  result  in  division  by  zero.  For  this  sit- 
uation, the  normalized  peak  value  of  the  modal  histor>  approxima- 
tion is  used  in  the  table.  The  invert  SSI  radial  maximum  stress  is 
an  example  of  this  situation. 

In  an  exactly  parallel  manner,  the  values  given  in  the  row 
labeled  MIN  are  computed  except  these  results  are  obtained  by  deal- 
ing with  the  minimum  values  of  the  SAMSON  and  modal  histories. 

The  first  row,  labeled  SRSS,  gives  the  results  from  performing 
square  root  of  sum  of  squares  (SRSS)  processing  on  the  data.  The 
last  column  gives  the  SRSS  of  the  nondimensional  crown  thrust.  The 
other  columns  are  determined  by  computing  the  SRSS  of  the  differ- 
ence between  the  modal  history  approximation  (for  the  number  of  modes 
indicated  in  the  column  heading)  and  the  SAMSON  response  history. 

The  desired  result  here  is  a ratio  of  0.0,  which  indicates  that 
there  is  no  difference  between  the  approximate  modal  history  curve 
and  the  SAMSON  time  history  curve. 

Obviously  the  desired  result  of  the  MODAL  HISTORY/ INPUT  VALUE 
ratios  is  1.0  for  the  MIN  and  MAX  rows  and  0.0  for  the  SRSS  row. 

A crown  thrust  history  based  on  modes  0 through  3 gives  a reasonably 
good  approximation  to  the  most  important  peak  values  of  the  response 
and  an  acceptable  SRSS  ratio.  Additional  modes  do  not  significantly 
change  the  peak  values  or  the  SRSS  ratio. 
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Some  of  the  ratios  show  some  seemingly  disparaging  results  such 
as  the  crown  SSI  radial  stress  MAT  row.  The  inclusion  of  modes  0 
through  6 for  the  approximation  to  this  variable  does  not  yield  the 
desired  result  of  approximately  1.0.  However,  an  examination  of 
figure  27c  reveals  the  fact  that  the  maximum  radial  stress  (tension) 
is  so  small  that  even  an  error  of  an  order  of  magnitude  is  not 
significant . 

4.  ANALYSIS  OF  CYLINDRICAL  SHELL  TO  SSI  LOADING 

This  analysis  follows  closely  that  presented  in  reference  9. 

It  presents  the  response  analysis  of  an  elastic  cylindrical  shell 
to  a transverse  shock  wave.  For  this  research  investigation  the 
shock  wave  can  be  represented  by  the  SSI  radial  and  shear  stresses 
such  as  those  computed  in  the  dynamic  SAMSON  solutions . The  effect 
of  viscous  damping  of  the  cylinder  has  been  added  for  this  analysis. 
The  dynamic  modal  response  analysis  is  employed  for  the  solution, 
and  both  extensional  and  inextensional  effects  are  considered.  The 
extensional  effects  are  frequently  neglected  for  analyses  of  this 
type;  however,  they  are  important  during  the  early  time  engulfment 
process . 

The  basic  objective  of  this  analysis  was  to  see  if  the  SSI  modal 
loads  computed  from  the  dynamic  SAMSON  solutions,  when  applied  to  a 
cylinder  in  vacuo,  would  produce  the  same  modal  internal  force  re- 
sultants (moment  and  thrust)  that  were  computed  from  the  SAMSON 
results.  The  theoretical  aspects  of  the  modal  analysis  are  summar- 
ized in  the  following. 

a.  Displacements  (see  figure  12) 


w(9,t)  = q0(t)  qn(t)  + qn(t) 

n=l  l 


cos  n9 


v K(t) 
= L -a— 

n=l  L n 


- dn  qn(t)  sin  n6 


where  the  bar  notation  ( ) denotes  the  extensional  effects  and  the 

comparable  unbarred  variable  the  inextensional  effects.  Thus,  a and 

^n 


i 


qn  are  the  time  varying  inextensional  and  extensional  modal  ampli- 
tudes, respectively.  Define  Nj = n2-l  , C = t2/(12R2)  then 


d -*k  + i 

dn  2n  ' 2n  C + 2 


2 2N2 

n +1  ZN1  N1  2 

~ir  — 2~  C + — 2 C 
' n^  n 


1/2 


dn  " 


Note  that  dg  = 0 and  d^  =1. 


Exact  Circular  Frequencies 
2 

0) 

n 

-2 
U) 

n 


where 


a = 


e = 


- 2m2 
e f n 

EC/(PR2) 


f = eC 

P = Y/S.  mass  density 

c.  Approximate  Frequencies 

If  n < 5 and  AR2  >>  I or  12  (I/t)2  » 1 

2 „ !sl  Nin2  , £(l-n2)2n2 

“n  ptR  (n2+i)  n2+l 

—2  ^ ^ 2 2 
“n  “ p|s  <"  +1>  * +1> 

d.  Radial  and  Tangential  Shear  Loads 


ar(0.  t)  = Y.  °rn<'t^  COS 


n9 


n=0 


(23a) 

(23b) 


= a - y a2  - b n > 1 

(24) 

= a + ~\j  ^ - b n > 0 

(25) 

0.5  e(n2+l)  +0.5  f N2 

(26) 
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xr9(0,t)  = ^ Tr0n sin  n9 


n=l 


e.  Equations  of  Motion 


(27) 


q‘  + 2 t to  q + to  q 
Hn  ^n  n Mn  n Hn 


arn(t)  ~T  Tr0n 
n 

m 

n 


(t) 


n > 1 (28) 


q + 2 £ to  q +oo  q ~~ 
Hn  n n Mn  n Mn 


arn(t)-  Tr0n(t) 
an 


m 


n 


where 


n > 0 (29) 


mn  " pt(1  +7Z) 

(30p) 

U 

n 

m = pt(i  L d2) 

(30b) 

n r 

^ ; dt 

?n  = Fraction  of  critical  damping 

Note  that  1/?  = 0 for  n=0 . 

n 

To  get  the  solution  for  the  separate  effects  of  the  radial  (R) 
and  tangential  (T)  shear  loads  the  following  equations  are  presented. 

Equations  of  Motion  for  Radial  Loads 
% + 2sn  “n  % + “n  % * " "i1  (31a> 

%+2=n;J5+;n^'  - ”rn(t)S„’  n>-°  (31b> 


Equations  of  Motion  for  Shear  Loads 
’n  + 2in  "n  + “n  ’n  ’ Tr9n(t)/dn/mn'  "i1  (31c) 

+ “n  ^ + “n  ’’n  ' Tr9n(c)/3n/V  "t1  (31d) 
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Note  that  for  equation  (31d)  the  n=0  mode  does  not  exist  since 
rrQn  = 0 and  1/3^  = 0 for  the  zero  mode.  The  complete  solution  is 
obviously 

R T 

qn  = qn  + qn  n >_  1 Inextensional  solution  (32a) 


qn  = 


n = 0 


q R + q T n > 1 

Mn  Mn  — 


Extensional  solution 


(32b) 


(32c) 


f.  Cylinder  Force  Resultants 


where 


M = — (w  + w") 
R 


T = — (w  - v' ) - — 


K = Ec1  = Ect3/12 

D = A E = E t 
c c 


( )•  = dLl 


(33a) 


(33b) 


(33c) 


(33d) 


The  radial  and  tangential  displacements,  equations  (21)  and  (22), 
can  also  be  written  as 


where 


where 


w(0 , t)  = w (t)  +/  w (t)  cos 
o i — ■ n 


wo(t)  = qo(t) 
Wn(t)  = qn(t)+qn(t) 


v(9 , t)  =^vn(t:)  sin  n0 


q (t) 

v (t)  = -2. - d q_  (t) 


(34a) 


(34b) 

(34c) 


(35a) 


(35b) 


Thus , 


w + w"  = wQ  + ^ (1-n^)  wn  cos  n6 
n=l 

00 

w-v'  = w^+  ) (w  -nv  ) cos  nQ 
o l — « n n 


From  equations  (33a)  and  (36) 


= 4 <w  + / (l-n^)w„  cos  nQ 

— / o n 


(38a) 


- - M%  + L (1‘n2)Wcos  110 

R L n=l 


M = M0+^Mn  cos  nQ  + 


Z M 

t—  n 


cos  nQ 


(38b) 


(38c) 


wnere 


M = Moments  from  extensional  effects 
n 

\=-^(l-n2)qn  n>0 

R 

M = Moments  from  inextensional  effects 
n 

y.  K ,,  2.  , 

Mn  ) qn  “i1 


(38d) 


(38e) 


Note  that  both  components  of  the  moment  are  zero  for  mode  1 since 
2 

1-n  =0.  From  equations  (37)  and  (33b)  obtain 


r oo 

T = — w + 4"  (w  -nv  ) cos  nQ  - 
5-  i o n n 

R - n=l  J 


(39a) 


n _ v-  _ nq 
- q + ) (q+q  - -r-2-  + nd  q ) cos  nQ 
B-  o n ^n  d„  n Mn' 

i\.  i n 

L n=l 


(39b) 
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T - T0+ZTn 


cos  n9  + 


Tt 

n 


cos  n9  - M/R 


(39c) 


The  last  term  (M/R)  is  usually  neglected;  however,  it  is  retained 
in  this  analysis. 


In  equation  (39c) 


T = Thrusts  from  extensional  effects 
n 


= | <ln<1+ndn)  + jp  (1-n2)  % 


= - q (1-n/H  ) + -4-  (1-n2)  qn  n>0  (39d) 

R n n RJ  n 


Note  that  1-n/d  = 1 and  2 for  n=0  and  1. 

n 


T = Thrusts  from  inextensional  effects 
n 


= | qn(l-n/dn)  + ^ (1-n2)  qn  n>  2 (39e) 


Also  note  that  d =n  for  n~l  thrust  T;  =0  and  this  has  been  used  to 
n x 


define  the  applicable  range  of  n for  equation  (39e) . 


The  nondimens ional  equations  of  internal  liner  forces  for  the 
separate  effects  of  radial  and  tangential  SSI  loads  can  be  written 
as  follows. 


g.  Nondimens ional  Moment  Modal  Amplitudes 


= Nondimens ional  modal  amplitude  for  SSI  °r 


pQR  load  only 


qR/p  , n=0  since  qR  = 0 for  n=0  (40a) 

n o ^n 


^ 2 J 0,  n=l  since  1-n  =0  for  n=l  (40b) 

* a (1-n) < 

R R R -R  ' 

n > 2 (40c) 

>po  Po  / 


= Nondimens ional  modal  amplitude  for  SSI  xrg 
load  only 


K / 1 2v 

- - -7=7 (1-n  ) 


0,  n<2  since  r^g^O  for  n=0 

and  l-ni. 2*0  for  n=l  (40d) 


T -T 

^n  , ^n  9 

— + — n > 2 

Po  Pq 


(40e) 


n _ 


Nondimens ional  Thrust  Modal  Amplitudes 
r T “R 


D , K n 2.  4n 
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RR  R RZ  Po 
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-R 
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(41a) 


(41b) 


B. ! + K 

^n;  R2R' 


-»  R 

.2,1  ^n 


n >2  (41c) 

(4ld) 

(41e) 


_5  11-— ' + 


/ -i  „ ^ \ n 1 D » , n * _i_  K,  ( i „ ^ \ nn 
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JPo 
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i.  Application  of  Method 

The  above  analysis  has  been  applied  to  two  of  the  SAMSON  solu- 

tions, cases  5 and  13.  A computer  code  was  written  which  used  the 
SSI  radial  and  shear  stress  modal  loads  that  were  computed  from  the 
SAMSON  solution  (see  section  III  for  derivation  of  these  modal 
histories) . These  modal  loads  were  used  in  the  equations  of 
motion,  equations  (31a),  (31b),  (31c)  and  (31d) , which  were  solved 
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numerically  by  a Runge-Kutta  technique  for  the  extensional  and  in- 
extensional  modal  displacement  histories  for  all  modes  up  to  n=7 . 
These  results  were  then  used  to  compute  the  nondimensional  modal 
thrust  and  moment  variations,  equations  (40)  and  (41).  Finally, 
the  liner  forces  computed  by  this  process,  which  are  referred  to 
as  the  analytical  liner  forces,  are  compared  to  those  computed  from 
SAMSON  results.  These  latter  liner  forces  are  referred  to  as  the 
SAMSON  liner  forces.  The  analytical  modal  thrusts  compared  fairly 
well  with  the  SAMSON  modal  thrust  for  both  cases  5 and  13.  The 
modal  moments  did  not  compare  very  well  for  some  of  the  modes, 
especially  mode  2. 

The  details  of  the  method  employed  to  make  the  comparisons  and 
to  measure  the  influence  of  the  separate  effects  of  the  radial  and 
shear  SSI  stresses  are  now  discussed.  Tables  VIII  and  IX  give 
the  vibrational  characteristics  of  modes  0 through  6 for  the  two 
cylinders  for  cases  5 and  13  (R/t  = 4 and  15) . The  outer  radius  and 
the  liner  thickness  used  for  these  cases  are 


R/t 

R,  inch  (m) 
t,  inch  (m) 


Case  5 


104.99  (2.667) 
26.25  (0.667) 


Case  13 
15 

84.37  (2.143) 
5.63  (0.143) 


These  values  are  consistent  with  the  2.0m  inside  tunnel  radius 
that  has  been  used  in  this  project.  In  the  tables  for  the  n=T 
mode,  a very  small  value  for  the  inextensional  frequency 

_ g 

(0.1  x 10  ) has  been  used  instead  of  0.0. 

The  column  headings  for  tables  VIII  and  IX  denote  the  follow- 


ing: 


Heading 
OMEGA- SQ. 
OMEGA- BAR- SQ 
FREQ. 

FREQ. -BAR 
T 


Parameter 


“n 

-2 

wn 

fn=o)n/2iT,  frequency  (Hz) 


n 


T =l/f  , period  (sec) 


Table  VIII 

MODAL  VIBRATION  PARAMETERS  (R/t = 4)  CASE 
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In  an  effort  to  efficiently  assess  the  effect  of  each  SSI 
modal  load  component  on  the  response  of  the  cylinder,  maximum  and 
minimum  analytical  values  of  the  moment  and  thrust  at  the  crown, 
springline  and  invert  were  computed  and  compared  with  the  corre- 
sponding peak  values  obtained  from  the  SAMSON  time  history  re- 
sults. These  comparisons  are  given  in  tables  X through  XIII 
and  are  similar  to  those  in  table  VII.  The  modal  forces  (thrusts 
and  moments)  used  for  table  VII  were  computed  directly  from  the 
SAMSON  tunnel  liner  force  histories  whereas  tables  X through  XIII 
were  computed  from  the  response  of  the  cylinder  in  vacuo  to  the  SSI 
(SAMSON)  modal  loads  for  tables  IX  through  XII.  Furthermore,  the 
time  of  occurrences  of  the  peak  moment  and  thrust  for  both  the 
analytical  and  the  SAMSON  results  are  given  in  tables  X through 
XIII  to  compare  this  important  aspect  of  the  problem. 

The  comparisons  given  in  these  tables  are  described  by  discus- 
sing the  "Crown  Thrust"  results.  First  consider  the  two  rows 
labeled  MAX-2  and  T-MAX-2.  The  last  column.  Exact  Response,  gives 
the  maximum  value  of  the  nondimensional  crown  thrust  (0.60422)  and 
the  time  of  occurrence  (0.0081  sec)  from  the  SAMSON  solution  for 
reference  point  1.  Knowing  the  analytical  modal  thrust  amplitude 
histories  (T^,  T^)  for  the  separate  effects  of  the  radial  and 
shear  SSI  loads,  the  crown  (Q^  = 3.75  deg)  thrust  history  is  calcu- 
lated for  a mode  by  mode  addition  of  these  modal  amplitude  histories. 
The  maximum  value  of  the  crown  thrust  was  computed  for  each  addi- 
tion of  a mode  to  the  analytical  history,  and  then  this  maximum 
thrust  is  normalized  by  dividing  by  the  SAMSON  maximum  value. 

These  are  the  normalized  values  that  are  presented  in  the  MAX-1 
and  MAX-2  rows  for  the  columns  labeled  MODE=0  through  6.  As  indi- 
cated in  the  table,  the  -2  row  includes  the  response  from  both  the 
radial  and  the  shear  SSI  modal  loads  for  modes  0 through  N where  N 


132 


TABLE  X 

COMPARISON  OF  PEAK  THRUSTS  FROM  ANALYTICAL  MODES  WITH  SAMSON  RESULTS,  CASE 
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CASE  XIII 

COMPARISON  OF  PEAK  MOMENTS  FROM  ANALYTICAL  MODES  WITH  SAMSON  RESULTS,  CASE  13 


JHIS  PAGE  IS  BEST  ftUALITT  PRACTICABLE 

® wr-M*  PUtfh  TO  HDC 


m 


MAX-1  0.000000  O.tioOOOO  .OJ1529  9.044965  11.145795  11.679975  11.579772 

-2  O.OOOOOO  O.uoPOOO  7 . 0525t>  1 12.38flUbi  11.317717  11. 342194  11.545565  ,43955t-02 

MAX-1  0.0000000  0.0000000  .0118000  .0294000  .0282000  .0294000  .0294000 

-2  O.000Q0O0  0.0O0OQ90  ,0294000  ,0294000  .0259000  .0294000  .0294000  .0173000 


denotes  the  column  label  mode  number.  The  -1  row  does  not  include 
the  SSI  shear  load  effects  for  the  current  mode  N.  For  example, 
using  modes  0 through  5 to  obtain  the  results  given  in  the  MODE  5 
column,  the  MAX-1  result  is  obtained  from 


Max  Tn  , = Max^  > 


f *• 


■0,5 


L~‘ 

'-n=0 


TR(t . ) + TT  (t . ) 
nv  iy  n iy 


cosn  9, 


™R 


+ (t^)  cos5  9^>,  i=l,L 


where  the  ( ) notation  indicates  the  nondimensional  magnitude,  n 
is  the  mode  number,  i is  the  discrete  time  step,  and  there  are 
L time  steps.  This  maximum  divided  by  the  SAMSON  maximum 
(0.60422)  would  give  the  result  0.987907.  The  T-MAX-1  and  -2 
rows  give  the  time  that  the  respective  analytical  maximums  occur. 


In  a similar  manner,  the  first  four  rows  (labeled  MIN-1, 
MIN-2,  T-MIN-1,  and  T-MIN-2)  are  obtained  except  the  minimum  values 
of  the  modal  histories  are  computed  and  compared  with  the  SAMSON 
solution  minimum  nondimensional  crown  thrust  (-0.040613)  and  time 
of  occurrence  (0.015225  sec). 


j . Discussion  of  Results 

It  was  anticipated  that  the  cylinder  moment  and  thrust  modal 
results  obtained  from  this  anaysis  would  be  close  to  those  obtained 
from  the  SAMSON  results;  however,  an  examination  will  show  that 
the  comparison  is  fair  for  the  thrusts  and  relatively  poor  for  the 
moments  for  case  5.  The  moment  comparison  is  especially  bad  for 
case  13.  Just  as  with  the  table  VII  results,  the  desired  ratio 
of  the  peak  maximum  and  mimimum  analytical  response  to  the  corre- 
sponding SAMSON  peak  is  1.0.  In  general,  the  minimum  thrust  ratios 
are  not  very  close  to  1.0;  however,  since  these  ratios  are  the  ten- 
sion thrust  which  are  relatively  small  except  at  the  invert  for 
case  5,  large  errors  are  not  significant.  That  is,  a large  error 
in  a small  magnitude  number  still  gives  a small  magnitude  number. 

It  also  appears  that  large  errors  exist  for  the  time  of  occurrence 
of  some  of  the  analytically  determined  peaks  compared  with  the  SAM- 
SON times;  however,  these  actually  represent  only  minor  errors. 
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For  example,  the  case  5 maximum  invert  thrust  for  the  SAMSON  solu- 
tion occurs  at  t=0.0096  sec,  whereas  the  maximum  occurs  at  approx- 
imately 0.026  sec  when  modes  0 through  2 (or  more)  are  included  in 
the  analytical  determination  of  the  invert  thrust.  However,  an  exam- 
ination of  figure  29a  shows  that  a second  large  compression  thrust 
occurs  at  approximately  0.025  sec  which  is  almost  equal  to  that  at 
t=0.0096  sec,  and  this  explains  this  seemingly  large  discrepancy. 

Even  though  the  comparisons  of  the  analytical  results  with 
the  SAMSON  solution  are  not  as  good  as  anticipated,  it  appears 
that  the  SSI  modal  loads  0 through  2 are  the  primary  contributors 
to  the  cylinder  thrust,  and  mode  2 is  the  significant  loading  caus- 
ing cylinder  moments.  As  previously  indicated,  the  desired  result 
was  to  have  the  analytically  determined  modal  thrust  and  moment 
histories  essentially  the  same  as  the  SAMSON  results.  A compari- 
son of  these  modal  histories  for  case  5 for  the  thrust  modes  0, 

1,  and  2 and  the  moment  modes  2 and  3 is  given  in  figures  30  and 
31.  The  thrust  modes  compared  fairly  well;  however,  relatively 
large  differences  exist  for  the  moment  modes  even  though  the  gen- 
eral time  variations  of  the  curves  are  similar,  i.e.,  the  analyti- 
cal moment  magnitudes  are  approximately  50  percent  higher  than  the 
SAMSON  magnitudes. 

The  thrust  and  moment  modal  comparisons  for  case  13  are  shown 
in  figures  32  and  33.  Just  as  with  case  5,  the  thrust  modes  com- 
pare fairly  well;  but  the  important  mode  2 moment  comparison  is 
very  poor,  and  this  is  the  major  error  which  creates  the  large  mo- 
ment ratios  in  table  XII I. 

A separate  ai  ysis  of  the  poor  comparison  of  the  analytical 
moments  with  the  SAMSON  results  was  undertaken  to  try  and  resolve 
this  problem.  This  study  was  not  completed,  due  to  fiscal  re- 
straints on  this  research  project;  however,  the  following  observa- 
tions are  made.  For  several  reasons  it  is  suggested  that  the 
SAMSON  determined  cylinder  thrusts  and  moments  are  more  reliable 
(or  accurate)  than  the  SSI  stresses.  Two  of  these  reasons  are: 

(1)  The  tunnel  liner  is  modeled  in  considerable  detail 
(four  element  layers  through  the  thickness)  com- 
pared to  the  soil  mesh  immediately  adjacent  to  the 
cylinder . 
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Figure  30a  Comparison  of  Analytical  Thrust  Mode  Variation 
with  SAMSON  Result,  Case  5 
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Figure  30b. 
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Comparison  of  Analytical  Thrust  Mode 
Variation  with  SAMSON  Result,  Case  5 
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Figure  31.  Comparison  of  Analytical  Moment  Mode  Variation 
with  SAMSON  Result,  Case  5 
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Figure  33.  Comparison  of  Analytical  Moment  Mode  Variation 

with  SAMSON  Result,  Case  13 


(2)  The  stresses  in  the  soil  elements  immediately  adja- 
cent to  the  cylinder  are  not  truly  the  SSI  stresses 
because  the  soil  mass  dynamic  effects  were  not  elimi- 
nated. That  is,  the  discrete  mass  at  the  SSI  nodes 
contain  contributions  from  both  the  soil  and  the 
cylinder  elements  that  are  attached  to  the  nodes.  For 
the  flexible  liner  (case  13)  the  soil  elements  are 
the  major  contributor  to  the  mass  at  the  SSI  node 
points.  The  stresses  in  the  soil  elements  immediately 
surrounding  the  structure  are  properly  used  to  find 
the  loads  applied  to  the  SSI  interface  nodes.  How- 
ever, it  is  clear  that  the  dynamic  force  effects  of 
the  soil  mass,  e.g.,  inertial  loads,  must  be  removed 
to  obtain  the  true  SSI  stresses.  Appendix  C gives 
a further  discussion  of  this  concept. 

Accepting  the  premise  that  the  SAMSON  thrusts  and  moments 

are  essentially  correct,  it  is  possible  to  solve  for  the  SSI  modal 

stresses  from  the  internal  liner  force  resultant  nodal  histories. 

This  is  done  by  first  solving  for  the  modal  deformation  histories 

(qn  and  q ) using  equations  (40)  and  (41) . Then,  by  application 

of  numerical  differentiation,  the  time  derivatives  of  q and  q 

n n 

can  be  obtained  and  equations  (31a)  through  (31c)  will  yield  the 
modal  histories  of  the  SSI  stresses. 

The  above  inverse  method  was  used  for  cases  5 and  13;  and 
even  though  the  results  are  preliminary,  they  imply  that  the 
mode  2 radial  and  shear  stresses  should  be  higher,  i.e.,  more 
positive.  A major  effect  occurs  on  the  case  13  moments  when  the 
higher  mode  2 radial  stresses  computed  from  tne  inverse  process 
are  used.  These  revised  radial  stresses  reduce  the  M2  amplitudes 

(make  them  more  negative)  which,  when  coupled  with  the  mode  2 
T 

moments  M2  from  the  SSI  shear  stresses,  yield  essentially  the 

same  total  mode  2 moment  variation  that  was  obtained  from  the 

SAMSON  results.  This  is  illustrated  in  figure  33  by  showing  the 

magnitudes  of  M2,  M£  at  t= 0.015  sec  as  computed  from  the  inverse 

process  (denoted  by  "inverse  method") . Also  shown  on  this  figure 
T R 

are  the  M2,  M2  magnitudes  at  t=0.015  sec  which  were  computed  from 
the  SAMSON  mode  2 SSI  radial  and  shear  stresses  (denoted  by  ORIG- 
INAL METHOD) . The  combination  of  these  two  moments  yields  a point 
on  the  "analytical"  curve.  Figure  34  shows  the  time  history  of 
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Figure  34.  Nondimensional  Mode  2 Moment  Variation  Computed  from  SAMSON  Mode  2 
SSI  Shear  and  Radian  Stress,  Case  13. 
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the  and  curves  as  computed  from  the  SAMSON  mode  2 SSI  radial 

and  shear  stresses.  Combining  these  two  curves  results  in  the 

dotted  curve  which  is  the  same  as  the  curve  denoted  "analytical" 

in  figure  33.  Figure  35  shows  the  nondimensional  mode  2 moment 

variation  that  was  computed  from  the  SSI  shear  and  radial  stresses 

T 

that  were  derived  by  using  the  inverse  method.  That  is,  the  M~ 

R ^ 

and  M2  curves  in  figure  35  are  based  on  the  SSI  mode  2 stress 

histories  calculated  from  the  SAMSON  mode  2 moment  and  thrust  his- 
tories as  discussed  above.  These  mode  2 SSI  stress  histories  are 
shown  in  figure  36  and  for  comparison,  the  SAMSON  SSI  stresses  are 
also  given.  As  previously  indicated,  the  SSI  stresses  determined 
from  the  inverse  method  are  more  positive  than  the  SAMSON  mode  2 
stresses  as  indicated  in  figure  36.  Referring  to  figure  35,  com- 
bining  the  M^  and  M^  (inverse  method)  moment  variations  yields 
essentially  the  original  mode  2 SAMSON  moment  variation  which  is 
the  desired  result.  In  a similar  manner,  the  mode  2 thrust  vari- 
ation computed  from  the  inverse  method  SSI  stresses  (figure  36) 
is  the  same  as  the  SAMSON  mode  2 thrust.  Similar  results  were  ob- 
tained for  the  other  inodes  when  the  inverse  method  was  used.  This 
indicates  that  the  computer  coding  associated  with  the  original 
analytical  method  is  correct  since  this  same  coding  was  used  to 
compute  modal  moments  and  thrusts  from  modal  SSI  shear  and  radial 
stresses  after  replacing  the  SAMSON  SSI  modal  stresses  with  those 
determined  from  the  inverse  method. 

In  summary,  it  is  concluded  that  for  the  results  available 
from  this  investigation,  the  SSI  stresses  should  be  computed  by 
the  inverse  method.  That  is,  the  modal  SSI  stress  variation  should 
be  calculated  from  the  SAMSON  modal  thrust  and  moment  histories. 
However,  further  research  studies  should  be  undertaken  to  verify 


this  conclusion. 


Figure  35.  Nondimensional  Mode  2 Moment  Variation 
Computed  from  Mode  2 SSI  Shear  and  Radial 
Stresses  Derived  by  Inverse  Method,  Case  13 


Figure  36a.  Nondimensional  Mode  2 SSI  Radial  Stress 
Variation  from  Inverse  Method,  Case  13 


SECTION  IV 


CONCLUSIONS  AND  RECOMMENDATIONS 

1.  CONCLUSIONS 

The  overall  objective  of  the  research  investigations  described 
in  this  report  was  to  develop  a simplified  method  for  computing 
loads  on  buried  structures.  The  results  discussed  in  Section  III 
indicated  that  the  soil  structure  interaction  (SSI)  loads,  as  de- 
termined from  the  dynamic  loading  finite  element  analyses,  may  not 
be  a sufficiently  reliable  data  base  upon  which  to  base  a new  or 
more  comprehensive  simplified  loading  technique.  Even  though  re- 
liable SSI  loads  are  not  available  some  meaningful  conclusions  can 
be  drawn  from  the  results  of  the  study.  The  conclusions  associated 
with  the  dynamic  loading  results  are  based  on  the  assumption  that 
the  detailed  results,  such  as  modal  time  history  variations,  may 
be  subject  to  question;  but  the  overall  effects  on  the  SSI  loads 
and  the  tunnel  liner  internal  forces,  from  varying  the  independent 
parameters  should  be  fairly  reliable.  For  example,  if  the  peak 
thrust  goes  down  with  an  increase  in  structure  flexibility,  the 
trend  is  essentially  correct  even  though  the  precise  reduction  may 
be  subject  to  question. 

The  major  conclusion  drawn  from  the  project  is  that  the  concept 
of  using  a few  of  the  SSI  modal  loads  to  arrive  at  a simplified 
method  of  loading  appears  to  have  a great  deal  of  merit.  In  fact, 
it  appears  that  it  is  generally  sufficient  to  use  only  modes  0 
through  2 to  obtain  a sufficiently  detailed  definition  of  the  load- 
ing to  satisfactorily  compute  the  internal  tunnel  liner  forces. 

In  some  cases,  it  may  be  necessary  to  include  mode  3 to  accurately 
obtain  the  bending  moments  around  the  tunnel  liner. 

Some  of  the  conclusions  that  can  be  drawn  from  the  static  and 
dynamic  investigations  undertaken  during  this  project  are  presented. 

In  this  discussion,  the  terms  "flexible"  and  "stiff"  refer  to  the 
tunnel  liner  with  radius  to  thickness  ratios  of  15  and  4,  respec- 
tively. Also,  "shallow"  and  "deep"  refer  to  depths  of  overburden 
above  the  tunnel  crown  of  0.25  and  2.0  tunnel  diameters,  respectively. 
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(1)  For  static  ground  surface  loads,  the  SSI  radial  and  shear 
stress  loads  acting  on  a buried  cylinder  are  not  very  sensitive 

to  the  depth  of  burial . 

(2)  For  dynamic  airblast  ground  surface  loads,  the  SSI  loads 
on  a buried  cylinder  are  significantly  reduced  as  the  depth  of 
burial  increases.  This  is  due  primarily  to  the  attenuation  with 
depth  of  the  peak  free  field  stress  for  the  nonlinear  soil  material 
models  (sand  and  cemented  sand)  considered  in  this  study.  This  at- 
tenuation with  depth  is  due  to  the  hardening  stress-strain  law 
employed.  Also,  a short  duration  airblast  load  (t^ = 0.015  sec) 
results  in  considerably  more  peak  stress  attenuation  than  the 

long  duration  loading  (^  = 0.06  sec)  considered  in  this  study. 

(3)  An  increased  failure  strength  of  the  soil  (i.e.,  cemented 
sand  versus  cohesionless  sand)  does  not  essentially  modify  the  SSI 
loads  on  a buried  cylinder  for  dynamic  airblast  loads. 

(4)  For  static  ground  surface  loads,  the  buried  cylinder  crown 
and  invert  thrusts  are  approximately  the  same,  and  the  springline 
thrust  is  two  or  three  times  larger  than  the  crown  and  invert 
thrust.  For  a shallow  depth  of  burial,  the  invert  thrust  is  slight- 
ly higher  (maximum  of  24  percent  for  the  stiff  tunnel)  than  the 
crown  thrust.  As  the  depth  of  burial  increases,  the  crown  thrusts 
increase  until  they  are  essentially  equal  to  the  invert  thrusts. 

This  occurs  at  approximately  one  diameter  of  overburden  for  the 
flexible  tunnel  and  two  diameters  for  the  stiff  tunnel. 

(5)  For  dynamic  airblast  loads,  the  crown  thrust  is  30  to  120 
percent  higher  than  the  invert  peak  thrust  for  the  stiff  tunnel; 
however,  for  the  flexible  tunnel,  the  crown  thrust  is  more  nearly 
equal  to  the  peak  invert  thrust.  For  the  shallow  buried  stiff  tun- 
nel, the  springline  peak  thrust  is  slightly  higher  (10  to  30  percent) 
than  the  crown  thrust;  otherwise,  the  springline  thrust  ranges  be- 
tween 50  and  150  percent  higher  than  the  crown  thrust.  The  peak 
thrusts  throughout  the  tunnel  liner  are  considerably  reduced  with 
depth  as  the  duration  of  loading  is  reduced  for  both  the  stiff  and 
the  flexible  tunnels. 
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(6)  For  static  ground  surface  loads,  the  crown  moment  is  some- 

i what  larger  (47  percent  higher)  than  the  peak  invert  moment  for 

the  most  flexible  shallow  buried  tunnel;  however,  the  crown  moment 
is  nearly  the  same  (7  percent  lower)  for  the  stiffest  shallow 
buried  tunnel  liner.  As  the  depth  of  burial  increases,  the  crown 
and  invert  moments  approach  equal  magnitudes.  This  result  occurs 
at  a depth  of  burial  of  one  tunnel  diameter.  The  invert  and  spring- 
line moments  are  essentially  of  equal  and  opposite  magnitudes  for 
all  depths  of  burial  and  tunnel  liner  flexibilities  considered. 

(7)  For  dynamic  airblast  loads,  the  largest  bending  moments 
occurred  at  the  invert  for  all  cases  considered.  This  rather  sur- 
prising result  corresponds  to  the  limited  experimental  results  avail- 
able at  this  time  (Ref.  58,59,60).  The  largest  nondimensional  mo- 
ments for  the  dynamic  load  cases  were  always  smaller  than  the  static 
solution  case  for  the  corresponding  R/t  ratio;  however,  the  invert 
peak  moments  approached  the  static  results  for  the  shallow  buried 
stiff  tunnel  liner  cases.  For  the  shallow  buried  stiff  tunnel, 

the  peak  dynamic  invert  moment  is  60  to  100  percent  higher  than 
the  crown  moment.  For  all  other  depth  of  burial  and  tunnel  stiff- 
ness combinations  considered  in  this  study,  the  peak  invert  moment 
was  20  to  50  percent  higher  than  the  crown  peak  value  except  for 
the  deep  buried,  flexible  liner  subjected  to  the  high  overpressure 
long  duration  loading.  For  this  case  the  invert  peak  bending  mo- 
ment was  80  percent  higher  than  the  crown  value.  In  most  cases, 
the  invert  peak  bending  moments  were  only  5 to  20  percent  higher 
than  the  springline  magnitudes  (but  of  opposite  sign) . In  some  of 
the  stiff  tunnel  cases,  the  magnitude  was  40  to  50  percent  higher. 

For  the  shallow  buried  stiff  tunnel  cases  subjected  to  short  dura- 
tion  airblast  loading,  moment  reversal  occurred  throughout  the  liner. 
The  maximum  reversal  (38  percent)  occurred  for  the  crown  bending 
moment  when  the  large  overpressure  short  duration  loading  was  used. 
That  is,  the  largest  negative  moment  (causing  tension  in  the  outer 
fibers  of  the  tunnel  liner)  was  38  percent  of  the  largest  positive 
moment  at  the  crown.  The  smallest  reversal  (16  percent)  for  the 
stiff  shallow  buried  cases,  occurred  for  this  airblast  load  case  also. 
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(8)  For  dynamic  airblast  loading,  an  increase  in  the  peak 
overpressure  has  a relatively  small  effect  on  the  peak  nondimen- 
sional  moments  and  thrusts  in  the  tunnel  liner.  The  largest  effect 
on  thrusts  is  an  increase  in  the  invert  compressive  nondimen- 
sional  thrusts  of  approximately  100  percent  when  the  peak  over- 
pressure is  tripled  for  the  shallow-buried  stiff  tunnel  with  the 
short  duration  loading.  All  other  changes  in  the  peak  nondimen- 
sional  thrust  with  an  increase  in  overpressure  are  relatively  small 
except  for  the  deep-buried  stiff  cylinder  subjected  to  the  long 
duration  loading  where  the  nondimensional  crown  and  invert  thrusts 
increase  by  approximately  50  percent.  The  largest  effect  on  bend- 
ing moments  is  a reduction  in  the  peak  nondimensional  moment  by 
approximately  40  percent  when  the  peak  overpressure  is  tripled  for 
the  shallow  buried  flexible  tunnel  liner. 

(9)  A fourfold  increase  in  the  duration  of  the  airblast  dy- 
namic loading  (from  15  to  60  msec)  has  the  effect  of  increasing  the 
peak  thrusts  and  moments  by  from  10  to  90  percent  for  the  shallow 
buried  tunnels.  For  the  deep  buried  tunnels,  the  increased  dura- 
tion of  loading  causes  large  increases  in  both  the  thrusts  and  the 
moments  which  generally  range  from  150  to  350  percent. 

(10)  Increasing  the  structural  flexibility  causes  a reduction 
in  the  crown  and  invert  nondimensional  thrusts  for  the  shallow 
buried  dynamic  load  cases;  however,  an  opposite  effect  is  observed 
for  comparable  static  load  cases.  The  nondimensional  bending 
moments  are  drastically  reduced  for  both  the  static  and  dynamic 
load  cases  as  the  tunnel  liner  flexibility  is  increased. 

2.  RECOMMENDATIONS 

This  research  project  produced  a great  number  of  numerical 
results  (36  cases)  from  dynamic  finite  element  solutions  which  mod- 
eled buried  cylinders  subjected  to  ground  shock  effects  of  airblast 
loading.  The  SSI  loads  on  the  cylinder  were  derived  from  the 
stresses  in  the  soil  elements  immediately  adjacent  to  the  cylinder. 
These  SSI  loads  were  then  used  in  a modal  analysis  to  compute  the 
response  of  a cylinder  in  vacuo,  and  the  internal  forces  (moments) 
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and  thrust)  computed  from  the  modal  analysis  were  compared  with 
those  determined  from  the  finite  element  solution.  In  general, 
the  thrust  compared  fairly  well,  but  the  moments  compared  poorly 
for  some  cases.  This  implies  that  either  the  SSI  stresses  or  the 
cylinder  forces  determined  from  the  finite  element  solution  are  in 
error.  At  this  point,  it  is  believed  that  the  SSI  stresses  are 
in  error  and  the  cylinder  forces  are  essentially  correct.  If 
this  is  true,  then  the  finite  element  determined  cylinder  forces 
can  be  used  to  calculate  the  SSI  stresses,  and  the  overall  objective 
of  determining  a simplified  method  for  defining  loads  on  buried 
structures  can  be  completed  by  using  the  results  from  the  36  cases. 

Briefly  stated,  the  simplified  loading  method  involves  the 
determination  of  the  important  modes  of  the  SSI  loads  and  an  analy- 
sis of  how  these  modal  histories  change  with  the  important  param- 
eters such  as  (1)  structural  stiffness,  (2)  depth  of  burial, 

(3)  soil  type,  and  (4)  airblast  waveform  peak  overpressure  and 
duration  of  loading.  It  appears  that  modes  0 through  2 provide 
the  most  significant  loads  on  the  structure;  however,  before  any 
simplified  loading  method  can  be  determined  from  these  SSI  modes, 
it  is  essential  to  have  confidence  in  the  SSI  modal  histories. 
Several  factors  have  been  identified  which  make  the  finite  element 
determined  SSI  loads  suspect  such  as  (1)  the  damping  stresses  are 
not  included  in  the  finite  element  stresses,  (2)  the  soil  mass  of 
the  elements  adjacent  to  the  structure  is  included  with  the  SSI 
nodes  and  (3)  a more  refined  mesh  may  be  required  to  accurately  de- 
termine the  SSI  loads  from  a finite  element  dynamic  analysis.  In 
summary,  it  is  recommended  that  this  research  study  be  continued 
and  completed  since  the  technique  of  defining  a simplied  loading 
procedure  by  using  a few  modal  loads  appears  to  have  great  merit. 
Specific  recommendations  are: 

(1)  Reliable  SSI  radial  and  shear  stress  modal  loadings  should 
be  determined  from  the  results  available  from  this  project.  Con- 
sequently, when  these  modal  load  histories  are  applied  to  the  cylin- 
drical tunnel  liner  structure  in  vacuo,  the  internal  forces  com- 
puted from  a dynamic  model  will  compare  favorably  with  those  ob- 
tained from  the  finite  element  (SAMSON)  solution. 


(2)  Once  the  SSI  modal  loadings  are  unquestionably  known  to 
be  without  error,  the  simplified  method  of  determining  dynamic 
loads  on  buried  structures  should  be  formalized.  This  effort  will 
require  the  determination  of  the  number  of  SSI  modal  loads  required 
and  the  time  history  details  (amplitude  and  duration  characteristics) 
of  each  mode  as  a function  of  the  structural  characteristics,  soil 
media,  depth  of  burial,  and  airblast  waveform.  The  time  details  of 
each  mode  can  be  simply  characterized  by  triangular,  piecewise 
linear  or  simple  sinusoidal  or  parabolic  waveforms  without  introduc- 
ing significant  errors  to  the  calculation  of  internal  tunnel  liner 
forces . 

For  example,  figures  37a  and  37b  show  the  mode  0 and  1 nondi- 
mensional  SSI  radial  stress  variations  for  cases  5 and  13.  These 
cases  are  identical  except  case  5 employs  a stiff  tunnel  lining 
(R/t=5)  and  case  13  a flexible  lining  (R/t=15).  The  other  important 
parameters  are:  peak  overpressure,  pQ=500  psi  (3.45  MPa);  airblast 
duration,  td=0.015  sec;  depth  of  burial,  h=0 . 5 R.  From  the  figures, 
it  is  seen  that  the  mode  0 radial  stress  can  be  simplified  by  a 
triangular  waveform  and  mode  1 by  three  linear  approximations  giving 
a large  amplitude  negative  triangle  followed  by  a smaller  amplitude 
positive  phase  triangle.  The  mode  0 triangular  approximation  is 
nearly  the  same  for  both  cases;  however,  the  peak  stress  for  the 
more  flexible  liner  (case  13)  is  slightly  higher  (15  percent)  and 
occurs  slightly  later.  Note  that  the  initial  slope  and  load  dura- 
tion of  the  mode  0 simplified  waveform  is  the  same  for  these  two 
cases.  For  the  mode  1 approximation,  the  initial  (negative)  peak 
stress  is  (65  percent)  higher  for  the  stiff  tunnel  (case  5)  than 
the  flexible  liner;  however,  the  positive  phase  peak  amplitude  can 
be  approximated  by  the  same  magnitude.  It  is  also  observed  that 
the  overall  load  duration  and  the  first  and  last  (third)  slopes  of 
the  mode  1 approximation  are  the  same  for  both  cases. 

Similar  simplified  SSI  stress  waveforms  can  be  determined  for 
higher  radial  stress  modes  and  the  SSI  shear  stress  modes  as  re- 
quired. It  has  been  indicated  that  modes  higher  than  mode  2 will 
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generally  not  be  required  to  accurately  determine  the  buried  struc- 
ture loads.  When  the  simplified  SSI  loads  are  characterized 
(as  a function  of  the  structure,  soil,  depth  of  burial  and  airblast 
loading  parameters)  the  reliability  of  these  buried  structure  load 
approximations  should  be  judged  by  computing  the  response  of  the 
cylinder  in  vacuo  to  these  loads  and  comparing  the  results  with 
those  obtained  from  the  original  unsimplified  modal  load  histories. 

(3)  The  results  of  the  dynamic  finite  element  solutions  con- 
sidered in  this  investigation  employed  two  different  nonlinear 
granular  type  soil  materials  with  a rather  high  yield  surface.  The 
results  from  the  solutions  using  these  two  materials  do  not  vary 
drastically,  and  this  implies  that  the  loads  on  a buried  tunnel 
liner  are  not  very  dependent  on  the  material  type.  However  a 
limited  number  of  results  were  obtained  during  the  project  using  an 
elastic  soil  media,  and  the  tunnel  liner  response  is  significantly 
different  than  with  the  nonlinear  media.  This  implies  that  the 
soil  media  characteristics  are  important  in  determining  the  loads 
on  buried  structures.  This,  of  course,  is  not  a surprising  con- 
clusion since  it  is  well  known  that  structures  buried  in  a competent 
material,  such  as  rock,  will  withstand  considerably  higher  ground 
surface  loading  than  when  buried  in  a poor  material. 

In  an  effort  to  quantify  the  relationship  between  SSI  loading 
and  soil  media,  it  is  recommended  that  solutions  be  obtained  using 
a clay-like  material  and  also  two  other  materials  which  would  simu- 
late very  good  and  poor  materials. 

(4)  It  was  observed  that  for  certain  of  the  dynamic  solutions, 
significant  tension  thrusts  occur  in  the  lower  portion  of  the  buried 
cylinder  liner  for  a portion  of  the  dynamic  solution.  In  some  cases 
these  tension  thrusts  exceed  the  peak  compression  thrust  that  occurs 
in  the  same  portion  of  the  liner  during  the  solution.  It  appears 
that  sizeable  tension  thrusts,  which  can  greatly  exceed  the  tensile 
capacity  of  a reinforced  tunnel  liner,  occur  for  shallow  buried 
stiff  cylinders  loaded  by  short  duration  ground  shock  waveforms. 
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The  soil  characteristics  are  also  an  important  aspect  of  this 
phenomenon,  since  for  the  stiff  tunnel,  only  minor  tension  thrusts 
occur  when  the  soil  is  considered  to  be  elastic.  Alternatively, 
large  invert  tension  thrusts  appeared  for  the  flexible  structure 
when  elastic  soil  was  used  instead  of  the  nonlinear  soil  media. 

The  tensile  thrusts  have  initially  been  attributed  to  the 
soil-structure  interaction  shear  stresses  that  occur  over  the  lower 
portion  of  the  cylinder  during  the  early  phases  of  the  loading. 

That  is,  the  initial  motion  of  the  cylinder  is  downward  in  re- 
sponse to  the  compression  wave  impinging  on  the  upper  portion  of 
the  cylinder.  This  motion  creates  SSI  shear  stresses  over  the 
lower  portion  of  the  cylinder  that  resist  the  downward  motion,  and 
these  shear  stresses  induce  sizeable  tensile  waves  in  the  invert 
portion  of  the  liner. 

Sizeable  tension  thrust  in  a buried  cylinder  has  not  been  pre- 
viously reported  by  other  analytical  or  experimental  studies  in- 
corporating realistic  soil  media  material  laws.  For  this  reason, 
it  is  recommended  that  this  phenomenon  be  investigated  in  more 
detail  so  that  the  mechanics  of  the  loading  are  understood  and  it 
will  be  known  what  parameters  (loading  waveforms,  structure  stiff- 
ness, soil  parameters,  etc.)  can  create  these  potentially  devas- 
tating structural  loads. 
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APPENDIX  A 


FREE  FIELD  VERTICAL  STRESS  WAVE  INVESTIGATIONS 

The  results  of  the  SAMSON  dynamic  finite  element  solutions  in- 
dicated that  when  the  cylindrical  tunnel  liner  was  buried  under 
two  tunnel  diameters  of  overburden,  relatively  drastic  reductions 
in  the  SSI  tunnel  loadings  occurred.  The  free  field  stress  wave 
study  presented  in  this  appendix  was  undertaken  to  quantify  the 
source  of  these  reduced  loadings  and  to  ensure  the  adequacy  of 
the  finite  element  model  used  to  represent  the  soil.  This  inves- 
tigation used  both  one-  and  two-dimensional  finite  element  meshes 
of  the  soil.  In  an  effort  to  see  whether  the  basic  soil  mesh 
used  in  the  buried  tunnel  liner  dynamic  solutions  was  creating 
some  of  the  attenuation  and  distortion  of  the  soil  stress  wave, 
the  two-dimensional  finite  element  mesh  was  made  up  by  filling  in 
the  tunnel  cavity  with  an  appropriate  mesh  of  soil  elements. 

Table  A1  gives  the  essential  parameters  for  the  14  free  field 
cases  discussed  herein.  The  results  of  this  study  ?^e  presented  in 
figures  A1  through  A9  and  show  that  the  primary  source  of  the  reduced 
SSI  liner  loadings  is  the  free  field  peak  stress  attenuation  with 
depth  for  the  two  nonlinear  soil  models  (materials  A and  B)  employed 
in  this  project.  That  is,  the  hardening  material  stress-strain 
variation  coupled  with  the  stiff  unloading-reloading  rate  creates 
rather  drastic  changes  in  the  airblast  pressure  waveform  as  it  prop- 
agates downward  into  the  soil  media.  The  most  pronounced  reduction 
in  the  peak  vertical  free  field  stress  with  depth  occurs  for  the 
short  duration  (0.015  sec)  surface  loading  case. 

The  finite  element  soil  meshes  were  made  up  as  follows.  First, 
the  representation  of  the  concrete  liner  was  removed  from  the 
two-dimensional  plane  strain  finite  element  model  and  the  resulting 
void  filled  in  by  126  soil  elements  as  illustrated  in  figure  A1 , 
thus  producing  a two-dimensional  soil  model  composed  of  803  rec- 
tangular and  triangular  soil  elements.  In  addition,  a simplied 
one-dimensional  model  was  constructed  which  consists  of  48  rec- 
tangular elements  representing  a vertical  column  of  soil. 


Table  A1 


FREE  FIELD  STUDY  RUN  MATRIX 


Case 

Mesh* 

P 

o 

(psi) 

Cd 

Soil 

Stress 

KL1 

-Strain  Load- 

Rates  (psi) 

KL2  ^L3 

Unload 

K 

u 

FF1 

5DSFF 

500 

0.015 

B 

16700 

38900 

166700 

166700 

2 

5DSFF 

1500 

0.060 

B 

16700 

38900 

166700 

166700 

3 

6 

500 

0.015 

B 

16700 

38900 

166700 

166700 

4 

6 

1500 

0.060 

B 

16700 

38900 

166700 

166700 

5 

6 

500 

0.015 

B-2  Load 

16700 

38900 

166700 

6 

6 

1500 

0.060 

B-2  Load 

16700 

38900 

- 

166700 

7 

6 

500 

0.015 

B-l  Load 

16700 

_ 

_ 

166700 

8 

6 

1500 

0.060 

B-l  Load 

16700 

- 

- 

166700 

9 

6 

500 

0.015 

_ 

16700 

38900 

- 

38900 

10 

6 

1500 

0.060 

16700 

38900 

- 

38900 

11 

6 

500 

0.015 

Elastic 

16700 

_ 

16700 

12 

6 

1500 

0.060 

Elastic 

16700 

- 

- 

16700 

13 

6 

500 

0.015 

A 

16700 

38900 

166700 

166700 

14 

6 

1500 

0.060 

A 

16700 

38900 

166700 

166700 

* 

Mesh 

5DSFF 

- Same 

as  Mesh  5DS  (h/R 

= 4)  with  cavity 

filled 

in 

with 

i soil 

media  elements 

Mesh 

6 

- One- 

dimensional  mesh  with  48 

elements 

1 psi  = 6.395  kPa 


Identical  load  cases  were  applied  to  both  of  the  models  and 
the  results  evaluated.  Results  from  the  two  models  compare  more 
closely  for  the  low  overpressure  (Pq=500  psi) , short  time  dura- 
tion (t^=0.015  sec)  loading  than  for  the  higher  overpressure 
(Pq=1500  psi),  long  time  duration  (^=0.60  sec)  case.  However, 
the  wave  shapes  and  peak  stress  values  obtained  using  the  detailed 
two-dimensional  mesh  for  the  soil  are  generally  in  good  agreement 
with  results  obtained  using  the  simplified  48  element  one-dimen- 
sional model  as  illustrated  in  figures  A2  and  A3.  This  indicates 
that  the  detailed  model  does  provide  a reliable  description  of  the 
soil  response  due  to  dynamic  loading. 
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The  soil  characteristics  provided  by  the  R0CKMT3  soil  model 
for  soil  B has  three  progressively  steeper  compressive  loading  rates 
and  an  unloading  rate  equal  to  the  steepest  loading  rate.  The  re- 
sponse of  the  soil  using  the  simplified  43  element  one-dimensional 
mesh  and  the  R0CKMT3  soil  description  is  shown  in  figure  A4 . For 
the  load  cases  considered  (PQ=500  psi  and  Pq=1500  psi)  the  result- 
ing stress  magnitudes  are  not  high  enough  to  require  use  of  the 
third  loading  rate.  Consequently,  if  only  the  first  two  loading 
rates  and  the  same  steep  unloading  rate  are  used  to  describe  the 
material,  the  soil  response  is  virtually  unaltered,  as  indicated 
by  figure  A5. 

By  using  the  first  two  loading  rates  and  an  unloading  rate 
equal  to  the  second  loading  rate  instead  of  the  steeper  third 
rate,  considerably  less  attenuation  of  the  stress  wave  occurs  as 
illustrated  in  figure  A7 . This  result  indicates  that  the  steeper 
unloading  rate  used  by  the  soil  A and  B models  is  primarily  respon- 
sible for  the  decay  of  the  propagating  stress  wave. 

Figure  A6  illustrates  the  result  of  using  only  the  first  load- 
ing rate  and  the  steep  unloading  rate.  For  the  case  of  low  over- 
pressure (Pq=500  psi)  only  the  first  load  rate  is  required;  there- 
fore, the  results  are  not  altered  from  the  PQ=500  psi  results  shown 
in  figures  A4  and  A5 . For  the  higher  overpressure,  the  second  load- 
ing rate  is  utilized;  and  as  expected,  a greater  attenuation  of  the 
stress  wave  is  realized  when  only  a single  loading  rate  is  used. 

Figure  A8  shows  the  one-dimensional  results  for  soil  A.  This 
material,  model  is  identical  to  soil  B except  the  yield  surface 
passes  through  the  origin,  i.e.,  there  is  no  cohesion  for  soil  A. 

The  results  for  soils  A and  B are  essentially  identical  as  seen 
by  comparing  figures  A4  and  A8 . 

Figure  A9  gives  the  results  for  an  elastic  soil  with  a bulk 
modulus  equal  to  that  used  for  the  first  loading  rate  of  soil  B. 

The  damping  parameters  for  this  material  are  also  the  same  as  those 
used  for  soil  B.  The  results  indicate  that  there  is  very  little 
attenuation  of  the  peak  stress  with  depth,  especially  for  the  high 
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overpressure  load  case  FF12.  There  is  also  a smaller  amount  of 
"shocking-up"  or  overshoot  of  the  peak  stress  at  shallow  depths 
(near  the  ground  surface)  in  comparison  to  the  soil  B solutions  as 
seen  by  comparing  figures  A4  and  A9 . 
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Figure  Al. 


Soil  Media  Mesh  to  Fill  In  Tunnel  Liner  Cavity 
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DEPTH/STRUCTURE  DIAMETER  (5.31  m) 


DEPTH/STRUCTURE  DIAMETER  (5.33  m) 


(e)  2D  Results  Case  FF1  (f)  ID  Results  Case  FF3 

Figure  A2  (concl) . Comparison  of  2D  and  ID  Free  Field  Vertical 
Stress  Waves,  Soil  B,  pQ  = 500  psi,  t^ = 0.015  sec 


DEPTH /STRUCTURE  DIAMETER  (5*33  ra) 


0 


(a)  2D  Results  Case  FF2  (b)  ID  Results  Case  FF4 

Figure  A3.  Comparison  of  2D  and  ID  Free  Field  Vertical  Stress  Waves 

Soil  B,  po=  1500  psi,  t(j  = 0.06  sec 
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DEPTH/STRUCTURE  DIAMETER  (5.33  m) 


DEPTH/ STRICTURE  DIAMETER  (5. 33  m) 


(f)  Maximum  Free  Field  Vertical  Stress,  2D  Results 

Figure  A3  (concl) . Comparison  of  2D  and  ID  Free  Field  Vertical 
Stress  Waves,  Soil  B,  pQ  = 1500  psi,  ^=0.06  sec 


Response  of  Bilinear  Loading  R0CKMT3  Soil 
= 16,700  KT  9 = 38,900  K = 166,700) 


Figure  A8.  Free  Field  Response  of  Trilinear  Loading  R0CKMT3  Soil  A 
(Kt  =16 , 700  , K,  2 = 38 ,900,  K, o=K  =166 , 700)  ID  Finrte  Element  Model 
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Figure  A9 . Response  of  Linear  Elastic  Soil  with  Damping 
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APPENDIX  B 

l 

LOADING  FUNCTIONS  ON  BURIED  CYLINDERS 
by  N.  M.  Newmark  Consulting  Engineering  Services 

As  an  airblast  wave  propagates  across  the  ground  surface,  it 
induces  a pressure  pulse  in  the  soil  which,  as  it  impinges  upon 
and  engulfs  an  underground  cylinder,  produces  loads  (pressures)  on 
the  cylinder.  Clearly,  as  this  soil  stress  wave  traverses  a buried 
cylinder,  the  pressures  produced  by  it  on  the  cylinder  vary  non- 
uniformly,  both  spatially  and  in  time,  around  the  cylinder. 

It  is  not  now  possible  to  define  these  soil-cylinder  inter- 
face pressures  with  a high  degree  of  confidence.  Qualitatively, 
they  should  be  influenced  by  a number  of  parameters  that  include, 
possibly  among  others,  the  depth  of  burial  of  the  cylinder,  the 
stiffness  of  the  cylinder  (in  both  flexure  and  in  radial  compres- 
sion) , the  properties  of  the  surrounding  soil  (including  local 
nonuniformities  in  it  around  the  cyliner) , the  period  (or  periods) 
of  vibration  of  the  cylinder  (or  cylinder-soil  system),  and  the 
size  (diameter)  of  the  cylinder  in  relation  to  length  of  the  free 
field  pressure  pulse.  Intuition  and  judgment  alone  provide  some 
insight  into  the  general  nature  of  the  influences  of  some  of  these 
parameters.  For  example,  as  the  depth  of  burial  increases,  the 
blast- induced  pressures  on  a buried  cylinder  should  be  reduced; 
as  the  stiffness  of  the  soil  increases  relative  to  that  of  the 
cylinder  which  it  surrounds,  a larger  proportion  of  pressure  that 
would  otherwise  be  resisted  by  the  cylinder  should  be  taken  by 
arch  action  in  the  soil  around  the  cylinder;  as  the  periods  of 
vibration  of  the  cylinder  become  longer,  the  nature  of  the  imposed 
dynamic  load  pulse  should  approach,  in  effect,  that  of  a statically 
applied  load. 

While  it  is  possible  to  discuss  these  probable  influences  in 
a qualitative  fashion,  they  are  not  yet  sufficiently  well  under- 
stood to  permit  one  to  define  them  quantitatively  with  a high  level 
of  confidence.  This  problem  is  presented  more  graphically 
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in  the  data  plots  (Ref.  1)  of  figures  B1  through  B8.  These  data 
were  taken  from  small  model  cylinder  tests  (generally  steel)  con- 
ducted under  controlled  conditions,  generally  when  buried  in  dense 
sand.  Inspection  of  these  data  plots  confirms  the  general  in- 
fluences of  some  of  the  parameters  discussed  above,  but  the  scatter 
of  the  data,  even  under  these  controlled  test  conditions,  is  so 
great  that  any  quantification  of  these  effects  is  almost  impossible, 
especially  insofar  as  flexure-producing  loads  are  concerned. 

Figures  Bl  through  B4  indicate  a general  reduction  in  circum- 
ferential thrust  with  increasing  depth.  However,  similar  plots, 
which  are  not  reproduced  herein,  fail  to  show  any  significant  re- 
lationship between  circumferential  thrust  and  cylinder  stiffness. 


Although  the  data  shown  in  figures  Bl  through  B4  are  quite 
scattered,  one  can  draw  a curve  as  shown  on  these  figures  that  at- 
tempts to  reflect  an  approximate  average  of  the  data  points.  The 
curve  shown  has  the  equation 


C 

a 


0.65  + 0.25 


<1.0 


= 0.65 


1.0 


in  which 

T^  = circumferential  thrust,  lb/inch 
p = free  field  pressure  at  depth  z = h , psi 
r = radius  of  cylinder,  inch 
hz  = depth  of  soil  cover  over  crown  of  cylinder 
D = 2r,  diameter  of  cylinder 


(B-l) 


Taken  from  "Design  of  Cylindrical  Reinforced  Concrete  Tunnel 
Liners  to  Resist  Air  Overpressure",  NCEL  Report  68.010,  prepared 
by  N.  M.  Newmark  Consulting  Engineering  Services,  June  1968. 
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Figure  B1 . Maximum  Thrust  Versus 


Figure  B2.  Maximum  Thrust  Versus  Depth  of  Burial  in  Sand,  p r'5/El=2-3 


Figure  B3 . Maximum  Thrust  Versus  Depth  of  Burial  in  Sand,  p rJ/El~3-4 


Crown  and  Invert  Moments  Versus  Normalized  Overpressure 

h =0  in  Sand 


Figure  B7 . Crown  and  Invert  Moments  Versus  Normalized  Overpressure 

h =D/4  in  Sand 


ao3ii 


Crown  and  Invert  Moments  Versus  Normalized  Overpressure 

h =D/2  in  Sand 


arching  factor.  This  equation  is  obviously  inadequate  since,  for 
a cylinder  of  any  given  diameter,  it  reflects  only  the  effects  of 
depth  of  cover  over  the  cylinder.  It  neglects  the  influences  of 
cylinder  stiffness  or,  perhaps  more  properly,  the  relative  stiff- 
nesses of  the  cylinder  and  the  surrounding  soil.  As  noted  earlier 
it  seems  extremely  unlikely  that  these  neglected  parameters  should, 
in  fact,  be  negligible. 

The  prediction  of  the  maximum  bending  moments  that  might  be 
experienced  by  a buried  cylinder  and,  hence,  the  prediction  of 
the  loads  that  would  produce  these  maximum  moments  is  fraught  with 
even  more  uncertainty  than  was  the  prediction  of  the  maximum 
thrusts.  These  greater  uncertainties  are  in  evidence  in  the  ex- 
tremely wide  scatter  of  the  experimental  data  shown  in  figures  B5 
through  B9 . Because  of  specialized  nature  of  these  data,  one 
should  be  extremely  careful  about  drawing  conclusions  from  them. 
Indeed,  they  may  show  little  more  than  judgment  would  tell  us  to 
expect;  namely,  that  for  flexible  cylinders  in  dense  materials 
only  modest  amounts  of  earth  cover  are  required  to  reduce  bending 
moments  to  relatively  negligible  amounts. 

Despite  the  scatter  of  these  data,  it  does  appear  that  the  mo- 
ment diminishes  quite  rapidly  as  the  cylinder  flexural  stiffness 
3 

(El/r  ) decreases  and  that,  for  the  range  of  data  points  available, 
the  moments  are  all  relatively  quite  small.  Furthermore,  although 
not  a clearly  defined  relationship,  there  is,  for  the  stiff er  cyl- 
inders, a trend  toward  reduction  in  moment  with  increasing  depth. 

As  suggested  by  the  above  observations,  one  of  the  primary  de- 
ficiencies of  the  data  contained  in  these  figures  is  the  absence 
of  data  for  flexurally  stiff  cylinders,  i.e.,  for  normalized  pres- 
sure values  less  than  1.0,  in  which  domain  it  is  noted  from  the 
plots  that  the  normalized  moments  appear  generally  to  be  increasing 
quite  rapidly. 
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These  data  were  obtained  either  from  static  tests  or  from  dy- 
namic tests  in  which  the  air  pressure  was  uniformly  applied  on 
the  ground  surface  above  the  cylinders.  Hence,  these  data  do  not 
properly  reflect  the  moments  that  would  be  induced  in  an  antisym- 
metrical  flexural  mode  during  passage  of  an  airblast  wave  across 
a cylinder,  thereby  loading  one  side  of  the  cylinder  before  the 
other  side  is  loaded.  There  is  reason  to  suspect,  especially  for 
very  shallow  buried  cases,  that  moments  associated  with  an  anti- 
symmetrical,  flexural  deformation  mode  (which  would,  or  might, 
result  from  the  horizontal  passage  of  an  airblast  wave  across  the 
ground  surface  above  a buried  cylinder)  may  be  greater  than  the 
moments  in  the  symmetrical  flexural  mode  (shortening  of  the  ver- 
tical diameter  and  lengthening  of  the  horizontal  diameter) . Be- 
cause of  the  manner  in  which  the  cylinders  were  tested,  the  data 
of  figure  B5  through  B8  contain  primarily  symmetric  flexural  modes. 

In  an  attempt  to  provide  a basis  for  a more  rational  repre- 
sentation of  the  influence  of  the  surrounding  soil  on  the  behavior 
of  a buried  cylinder,  studies  were  made  of  a half-cylinder  (i.e., 
an  arch)  in  which  the  soil  was  represented  as  a series  of  radial 
springs  and  added  masses  which  would  respond  compatibly  with  the 
cylinder  (or  arch)  as  it  deforms  under  load.  Unfortunately,  these 
studies  were  unsuccessful.  Results  obtained  therefrom  were  in 
substantial  disagreement  with  the  known  dominance  of  the  uniform 
compression  mode  of  response  for  fully  buried  flexible  arches  and, 
as  a consequence,  must  be  disregarded. 

However,  on  the  basis  of  these  studies,  it  appears  that  the 
soil  effects  on  the  loading  and  response  of  an  underground  arch 
(or  cylinder)  cannot  be  represented  simply,  at  least  not  as  direct 
modifications  of  the  structural  arch  mass  and  stiffness.  Indeed, 
additional  studies  suggest  that  a more  reasonable  approach  may  be 
to  model  the  soil  as  a soil  arch  which  surrounds  the  structural 
arch  and  whose  deformations  must  be  compatible  with  those  of  the 
structural  arch.  However,  to  define  such  a structural  arch-soil 
arch  model  rationally  and  simply  requires  major  research  that  is 
substantially  beyond  the  scope  of  this  effort. 


It  is  clear  that  the  present  inadequate  state  of  knowledge 
will  not  permit  analysis  of  a buried  cylinder  comprehensively, 
taking  into  account  properly  and  rationally  the  interaction  of  the 
cylinder  with  the  surrounding  soil.  Consequently,  it  seems  most 
reasonable  to  adopt  a procedure  which  has  at  least  the  virtue  of 
simplicity  and  which  yields  results  that,  in  the  limit,  appear  to 
be  reasonable  and  which  are  in  general  agreement  with  the  limited 
amount  of  test  data  that  are  available. 


Consistent  with  this  simplistic  approach,  an  underground  cylin- 
der can  be  assumed  to  respond  to  two  primary  modal  loadings, 
namely,  a uniform  compression  mode  and  an  anti symmetrical  flexural 
mode.  For  the  compression  mode,  the  pressure  on  the  cylinder  can 
be  taken  to  be  uniform  and  equal  in  intensity  to  the  free  field 
vertical  pressure  pgo  at  depth  hav>  where  hav  is  approximately 
(hz  + 0.2r),  except  that  the  rise  time,  t , is  increased  by  one- 
half  the  transit  time  of  the  blast  wave  across  the  cylinder.  This 

transit  time  is  equal  to  L /U,  where  L = D,  diameter  of  the  cylin- 

s s 

der  and  U is  the  shock  wave  velocity.  The  resulting  function  is 
shown  in  figure  B9(a).  These  pressures  are  reduced  to  account  for 
the  arching  effect  of  the  soil  by  the  factor  C , which  may  be 
estimated  from  equation  (B-l) . 


For  the  flexural  mode,  the  pressure  distribution  is  assumed 
to  be  antisymmetrical , with  uniform  inward  pressures  on  the  upper 
windward  and  lower  leeward  quarters  of  the  cylinder,  and  equal 
uniform  outward  pressures  on  the  other  two  quarter  segments.  The 
pressure  intensity  on  each  of  the  four  quarter  segments  is  the 
same,  and  is  assumed  to  vary  in  time  as  shown  in  figure  B9(b).  The 
factor  Cfz>  which  appears  in  this  pressure-time  function,  is  intro- 
duced to  provide  for  an  orderly  reduction  in  this  flexural  mode 
loading  as  soil  resistance  to  structural  deformation  in  this  mode 
increases  with  increasing  depth  of  cover.  Somewhat  arbitrarily, 
this  factor  is  taken  as 


Cfz  1,2 


'h  ' 
z 

r i 


(B-2) 
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which  is  never  taken  greater  than  1.0  nor  less  than  zero.  As  is 
evident  from  this  expression,  the  flexural  mode  is  not  considered 
to  exist  if  hz/r  is  greater  than  0.3.  This  is  clearly  an  arbitrary 
judgment  which  fails  to  acknowledge  that  the  resistance  offered 
by  the  soil  to  structural  deformation  in  this  mode  must  be  a func- 
tion not  only  of  the  depth  of  soil  cover  but  also  of  the  soil 
type,  density,  and  stiffness,  or  resistance  to  deformation. 

As  noted  earlier,  these  suggested  modal  loadings  are  obvious 
simplifications  of  a very  complex  problem.  Much  more  research 
involving,  perhaps,  finite-element  solutions  of  the  highly  non- 
linear and  inelastic  soil- structure  continuum  is  needed  to  assess 
the  reasonableness  of  this  approximation.  Desirably , ' in  view  of 
the  known  difficulties  involved  in  modeling  such  a system,  the 
results  of  such  theoretical  studies  should  be  checked  experimentally. 

A few  final  comments  in  this  regard  seem  appropriate : 

(1)  The  conditions  of  interest  are  those  that  correspond 
to  failure  of  the  cylinder,  recognizing  that  failure 
is,  itself,  difficult  to  define  in  many  cases.  But 
almost  certainly,  we  are  interested  in  the  behavior 
of  the  cylinder,  and  probably  of  the  surrounding 
soil,  well  beyond  the  range  of  linear  elastic  response. 

(2)  Unless  the  cylinder  is  quite  long  relative  to  its 
diameter,  the  effects  of  the  ends  of  the  cylinder  on 
its  behavior  under  applied  external  loads  should  be 
considered.  Limited  studies  have  indicated  that  the 
shell  effects  of  closed  cylinder  ends  can  increase 
markedly  the  ability  of  a relatively  short  cylinder 
to  withstand  flexural  loading. 

1.  UNDERGROUND  PRESSURE-TIME  FUNCTIONS 

a.  General  Philosophy  and  Concepts 

There  is  now  no  rational  basis  on  which  to  define  the  detailed 
character  of  the  air-induced,  underground  pressure  pulse  shape. 

It  is  clear  that,  as  a given  surface  overpressure  pulse  is  propa- 
gated vertically  downward  through  a soil  medium,  its  shape  will 
be  changed;  both  theoretical  wave  propagation  studies  in  nonelastic 
materials  and  field  test  data  confirm  this.  It  is  not  yet 
possible  to  define  the  highly  variable  in-situ  soil  dynamic 
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properties  within  the  limits  needed,  nor  is  it  possible  to  repre- 
sent the  wide  variety  and  the  infinite  possible  combinations  of 
soil  types  and  conditions  with  a mathematical  model  that  would 
permit  the  development  of  a theoretically  rigorous  underground 
pressure  pulse  shape. 

Similarly,  while  the  nuclear  field  tests  have  provided  useful 
information  in  regard  to  these  pulse  shapes,  the  data  scatter  to 
such  an  extent,  even  in  repeated  tests  in  the  same  soil,  and  are 
available  in  only  a few  relatively  atypical  soils,  that  detailed 
underground  pulse  shapes  cannot  be  predicted  with  confidence  from 
this  source  either. 

Nevertheless,  pressure  pulse  shapes  must  be  known  (or  estimated) 
in  order  to  investigate  the  vulnerability  of  buried  structures  to 
ground  surface  airblast  loadings;  hence,  it  is  suggested  that  they 
be  approximated  on  the  basis  of  the  following  criteria  which  are 
defensible,  at  least  in  general  terms: 

(1)  The  soil  is  assumed  to  be  a uniform  medium  of  in- 
finite depth.  No  attempt  is  made  here  to  account 
for  the  effects  of  a layered  system. 

(2)  The  total  impulse  that  is  applied  to  the  ground  sur- 
face must  be  propagated  without  diminution  through 
the  underlying  soil  medium.  Hence,  the  iso-impulse 
curves  (discussed  subsequently),  which  were  devel- 
oped using  a pseudostatic  computational  technique 

on  the  basis  of  impulse  conservation,  can  be  used  to 
identify  the  total  impulse  that  should  exist  under 
the  vertical  pressure-time  function  at  any  depth  for 
any  given  weapon  yield  and  surface  overpressure  level 
or  range. 

(3)  The  peak  vertical  pressure  must  attenuate  with  depth 
of  penetration.  The  mechanism  by  which  this  pressure 
attenuation  takes  place  is  not  fully  understood,  al- 
though it  clearly  must  be  influenced  by  the  shape 

of  the  surface  overpressure  pulse,  the  level  and 
distribution  of  pressures  over  the  ground  surface, 
and  probably  the  soil  properties.  The  first  two 
effects  are  represented  by  the  a factor  of 
equation  (B-3).  The  latter  effect  is  not  yet  well- 
defined,  but  evidence  indicates  that  it  will  gener- 
ally be  relatively  small  in  comparison  to  the  other 
two  effects;  hence  no  effort  was  made  to  include 
it  herein. 


where  z = depth  below  ground  surface,  and 


L - 2300  (Wmt)1/3  ps;0-5  <5§j) 

in  which  C,  in  feet  per  second,  is  the  seismic  velocity 
of  the  soil  and  is  the  nuclear  weapon  yield  in  MT. 

(4)  As  the  pressure  pulse  penetrates  into  the  soil,  its 
rise  time  to  peak  pressure  increases,  but  neither 
the  duration  of  this  rise  time  nor  the  shape  of  the 
pulse  during  this  period  of  time  can  be  predicted 
closely.  It  is  suggested,  in  the  absence  of  better 
information,  that  the  rise  time  to  peak  pressure 

at  a given  depth  be  taken  as  1 msec  plus  half  the 
transit  time  of  the  peak  pressure  to  that  depth. 

Rather  arbitrarily,  the  shape  of  the  pulse  during 
its  rise  is  taken  to  be  parabolic.  This  is  done 
in  an  attempt  to  recognize  the  extreme  improbability 
of  the  occurrence  of  a linear  rise  to  maximum  pres- 
sure. Available  acceleration  test  data  support 
the  argument  that  during  rise  the  shape  of  the 
velocity  pulse  and,  hence,  of  the  pressure  pulse 
is  in  fact  nonlinear. 

(5)  There  being  virtually  no  information  available  in 
regard  to  the  shape  of  the  pressure  pulse  beyond 
maximum  pressure,  it  seems  reasonable  to  assume 
that  the  total  duration  of  the  pulse,  excluding 
the  rise  time  just  discussed,  be  taken  as  nearly  as 
possible  to  be  equal  to  the  positive  phase  duration 
of  the  surface  overpressure  pulse.  Hence,  the 
duration  of  the  pulse  beyond  maximum  pressure  may 
be  taken  equal  to  the  surface  overpressure  duration 
except  when  forced  to  vary  from  this  in  order  to 
comply  with  the  previously  stated  criteria  of  im- 
pulse conservation,  peak  pressure  attenuation,  and 
rise  time  pulse  shape. 

b.  Overpressure-Induced  Vertical  Impulse 

(1)  Positive  Phase  Overpressure  Impulse  on  Ground  Surface 

The  overpressure  impulse  produced  on  the  ground  surface  is 
given  approximately,  but  adequately,  by  the  following  equation. 

I = 1.739  p psi-sec,  for  0 < p < 1000  psi 

= 4.125  psQ^ ' psi-sec,  for  pgQ  > 1000  psi  (B-4) 
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(2)  Vertical  Free  Field  Impulse  Below  Ground  Surface 

For  a surface  overpressure  impulse  of  200  psi-sec  or  greater 
(for  W=1  MT)  , the  iso-impulse  curve  is  approximated  closely  by 
a circle,  as  shown  in  figure  BIO,  whose  center  is  at  ground  zero 
and  whose  radius  is  equal  to  the  ground  range,  R , determined  to 
correspond  to  the  impulse,  I,  of  interest. 
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Figure  BIO.  Representation  of  Underground  Iso-Impulse  Curves 
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For  a surface  overpressure  impulse  equal  to  or  less  than  200 
psi-sec  (for  W=1  MT)  the  iso-impulse  curve  is  adequately  repre- 
sented, as  shown  in  figure  BIO,  by  a vertical  straight  line  seg- 
ment from  the  ground  surface  and  a circular  arc.  The  straight 
line  extends  from  the  ground  surface  at  a range  Rn  (corresponding 
to  the  surface  range  as  determined  for  the  impulse  of  interest) 
to  a depth,  z , that  is  given  by 

zn  = 200  I-^^,  for  5 _<  I _<  200  psi-sec  for  W=1  MT  (B-5) 


The  lower  part  of  this  iso-impulse  curve  is  represented  closely 
by  a circle  whose  center  is  directly  below  ground  zero  (R = 0)  at 
a depth  of  z , and  whose  radius  is  equal  to  R for  the  impulse  of 
interest.  Thus,  the  circular  arc  is  tangent  to  the  straight  line 
segment  at  a point  whose  coordinates  are  R^,  Zn* 

These  approximate  iso-impulse  curves  are  plotted  in  figures 
Bll  and  B12. 

(3)  Notation  Used  for  Underground  Pulse  Shape  Determinations 

The  following  notation  is  used  in  the  description  of  the  under- 
ground pressure  pulse  shapes: 

Pso  = peak  surface  overpressure,  psi 

p = ap  = attenuated  peak  vertical  pressure  at  depth  z,  psi 
where  a is  given  by  equation  (B-3) 

t = surface  overpressure  positive  phase  duration,  sec 

t'  = positive  phase  duration  of  a surface  overpressure  pulse 
having  a peak  value  of  p^,  sec 

t = travel  time  of  pulse  from  ground  surface  to  depth  z,  sec 

t = rise  time  to  peak  value  of  pressure  pulse  at  depth  z 
(taken  here  as  [0.001+  (1/2) t ],  but  no  greater  than 
(l/2)ts)  C 

I = positive  overpressure  impulse  on  ground  surface,  psi-sec 

I'  = positive  impulse  under  a surface  overpressure  pulse 
of  peak  value  pz>  psi-sec 
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Figure  Bll . Approximate  Underground  Iso-Impulse  Curves 

(High  Impulse  Region) 
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Figure  B12.  Approximate  Underground  Iso-Impulse  Curves 

(Lok  Impulse  Region) 
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= impulse  under  rise  time  portion  of  underground  pressure 
pulse,  psi-sec 

Iz  = positive  phase  impulse  at  depth  z,  psi-sec 

c.  Definition  of  Underground  Pulse  Shapes 

For  individual  or  hand-computational  purposes,  an  adequate 
underground  pulse  could  be  approximated  or  sketched  for  a particu- 
lar case  on  the  basis  of  the  criteria  stated  above,  but  such  a 
procedure  is  not  satisfactory  for  a general  computerized  program 
of  analysis.  For  a computer  analysis,  it  is  essential  that  the 
underground  pulse  shapes  be  so  described  that,  for  a give  i set  of 
conditions,  it  is  uniquely  defined.  Consequently,  the  fc  lowing 
procedure  was  developed  for  computer  solution  to  define  under- 
ground pulse  shapes  consistent  with  the  general  criteria  stated 
above  and  which  are  uniquely  defined  by  a specific  set  of  given 
conditions . 

For  a given  surface  overpressure  level,  weapon  yield,  depth 
below  ground  surface,  and  soil  profile,  proceed  as  follows: 

1.  Assume  a parabolic  rise  in  the  underground  pressure 
from  zero  to  a maximum  pressure  of  p = ap  in  a 

time  t . Then  I = 2/3  pz  t . (B-6) 

2.  Assume  the  pulse  to  decay  from  its  peak  pressure  to 
zero  in  a manner  defined  by  a surface  overpressure 
pulse  having  a peak  pressure  equal  to  p , a corre- 
sponding duration  of  t^  and  an  impulse  of  I'. 

Steps  1 and  2 yield  a pulse  as  shown  in  figure  B13(a)  . 

This  pulse  is  satisfactory  when  (Ir+I')  = I . 

3.  If  (Ir  + I*)  <I?,  add  a parabola  as  shown  in  figure 
B13(b;  to  the  decay  portion  of  the  pulse  to  yield 
the  pulse  shown  in  figure  Bl3(c),  for  which  the 
total  impulse  is  I . 

4.  If  (Ir  + I^)  > Iz , proceed  as  follows: 

(a)  For  I <_  I , use  a parabolic  rise  time  portion 
as  defined  in  step  1,  which  is  followed  by  a 
decay  from  peak  pressure  as  defined  in  step  2 
except  that  the  decay  duration  is  St’s  instead 
of  tg,  in  which  g is  less  than  1.0  and  is  given 
by 
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Xz  ~ Ir 


(B-7) 


This  yields  the  pulse  illustrated  in  figure  B13 (d) . 

(b)  For  Ij.  > 1/2  Iz , arbitrarily  assume  a rise  time 
t^.  = 1/2  tr,  and  use  a straight  line  increase 
to  p Compute  the  corresponding  rise  time 
impulse  to  be 


I'  = 
r 


I Pz  Cr 


(B-8) 


For  step  4(b): 


(1)  If  (I'  +1')  < I2  assume  the  decay  portion  of  the  curve 
to  berdefine?  by  steps  2 and  3 , except  that  the  para- 
bolic segment  of  step  3 has  an  area  of  [Iz  - (I^  + I')] 
instead  of  [Iz  - (Ir  + I')].  This  produces  the  pulse 
illustrated  in  figure  B13(e)  . 


(2)  If  (I^  + I')  > Iz,  but  1^  < 1/2  Iz,  use  a linear  rise 
time  portion  of  duration  t^.  as  above,  followed  by  a 
decay  from  peak  pressure  to  zero  in  a time  B'ts>  in 
which 

I - I' 

6'  = -Sjt  - ( 


6'  = 


(B-9) 


This  yields  the  pulse  depicted  in  figure  B13(f). 


(3)  If  (1^  + I')  > Iz , but  1^.  >_  1/2  Iz , assume  the  pulse 

to  be  triangular  in  shape , as  illustrated  in  figure  B13 (g) 

The  detailed  character  of  the  outlined  procedure  should  not 
be  interpreted  as  implying  that  the  resulting  pulse  shapes  are 
precise  in  the  sense  that  they  represent  the  correct  pulse  shapes. 
As  stated  earlier,  correct  pulse  shapes  cannot  now  be  predicted 
with  a high  degree  of  confidence  by  any  method.  The  procedure  pre- 
sented here  is  considered  to  yield  pulse  shapes  that  are  as  cor- 
rect as  can  be  developed  conveniently  for  computer  solution  on 
the  basis  of  currently  available  data.  The  rather  complex  charac- 
ter of  this  procedure  results  solely  from  the  necessity  of  spe- 
cifying uniquely  defined  pulses  for  computer  solution. 
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APPENDIX  C 


DYNAMIC  INTERACTION  BETWEEN  ELASTIC  MEDIA 

In  this  appendix  a simplified  problem  of  dynamic  interaction 
between  two  elastic  media  is  considered.  The  objective  of  this 
analysis  is  to  determine  the  relative  influence  of  the  soil  and 
the  structure  discrete  mass  distribution  on  the  SSI  forces.  For 
the  dynamic  finite  element  solutions  described  in  Section  III, 
the  stresses  in  the  soil  elements  adjoining  the  tunnel  liner 
structure  were  used  as  the  SSI  loads.  The  analysis  given  herein 
addresses  the  adequacy  of  this  approach. 

To  study  the  problem,  the  structure  is  modeled  by  a single- 
degree of  freedom  system  with  mass  m^  and  stiffness  k ^ . The  soil 
element  adjoining  the  structure  is  modeled  by  two  lumped  masses 
of  m2  each,  and  an  element  stiffness  of  k2  joining  the  two  mass 
points  as  shown  in  figure  Cl(a).  The  combined  soil  structure  system 
is  shown  in  figure  Cl(b).  The  undamped  equations  of  motion  of  the 
combined  system  is  given  by 


m-^+n^ 


k1+k2 


(C-l) 


where  u-^  and  ^ are  the  translational  degrees  of  freedom  as  shown 
in  figure  Cl,  superscript  dot  represents  time  derivative, 
and  F is  the  force  applied  to  the  mass  nu,.  To  simplify  the  analy- 
sis presented  here,  it  is  assumed  that  F is  applied  suddenly  and 
remains  constant  with  time. 

The  two  combined  frequencies  of  the  simplified  system  are  given 
by  the  characteristic  equation 


2 If 

kf  + k2  - oj  (m^  + n^)  k. 


- k^  = 0 


(C-2) 


For  each  mode  the  normalized  mode  shapes  are 


(k2  - ) /D 


k2/D 


(C-3) 


l 


r ' n 

AFWL-TR-78-6 

i 


(a)  FREE  BODY  DIAGRAM  OF  THE  IDEALIZED  STRUCTURE 
AND  SOIL  ELEMENT 


Figure  Cl.  Simplified  Soil  Structure  Finite  Element  Model. 
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where 

D = ^(111^2)  (k2-m2u)2)2  + m2k2J1/2 
The  response  of  the  system  is  given  by 


where 


p = (Fk2/D) (l-coswt) 


(C-4) 


(C-5) 


The  force  that  is  applied  to  the  structure  (F^)  is  the  force  that 
mass  m2  transmits  to  mass  m-^  at  the  soil-structure  interface.  This 
force  is  given  by 

= F2  - m2u^  (C-6) 

where  F2  is  the  force  in  the  soil  element  (k2) 

F2  = k^u^up 

Equations  (C-3)  to  (C-7)  yield 

2 2 

F^  l-(2-a)  /o)2)coso)t 
F^  1- coswt 

2 

in  which  oj2  = k2/m2 

If  the  applied  force  is  a periodic  function  like  F cosftt,  the 
above  ratio  becomes 


(C-7) 


(C-8) 


2 2 1 

hA-A 

cosftt  - 

r 2 •) 

2 - (JT") 

Od  U)2 

Co 2 

cosoit 


cosfit  - COSCdt 


(C-9) 


It  is  apparent  from  equations  (C-8)  and  (C-9)  that  F^  and  F2 
can  be  substantially  different.  Therefore,  in  general  it  is  not 
desirable  to  use  the  soil  stresses  in  the  elements  adjacent  to 
the  structure  (F2)  as  the  SSI  stresses  (F^). 
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5301  Bolsa  Ave.,  Huntington  Beach,  CA  92647 
ATTN:  Robert  W.  Halprin 

Merritt  Cases,  Inc. 

P.  0.  Box  1206,  Redlands,  CA  92373 
ATTN:  Tech  Lib 
ATTN:  J.  L.  Merritt 

Newmark,  Nathan  M. , Consulting  Engineering  Services 

1211  Civil  Engineering  Bldg,  Rm  B106A,  Univ.  of  Illinois 

Urbana,  IL  61801 

ATTN:  Nathan  M.  Newmark 

ATTN:  W.  J.  Hall 

ATTN:  J.  Haltiwanger 

Pacifica  Technology 

P.  0.  Box  148,  Del  Mar,  CA  92014 

ATTN:  R.  Bjork 

ATTN:  G.  Kent 


219 


Dept,  of  Defense  Contractors,  Cont'd. 


No.  of 
Copies 


1 

1 


1 

1 

1 

1 

1 


1 

1 


1 


1 


1 

1 


1 

1 


1 

1 

1 


1 

1 


Physics  International  Co. 

2700  Merced  St.,  San  Leandro,  CA  94577 
ATTN:  Doc  Con  for  Tech  Lib 
ATTN:  Doc  Con  for  Coye  Vincent 

R & D Associates 

P.  0.  Box  9695,  Marine  Del  Rey,  CA  90291 

ATTN:  Jerry  Carpenter 

ATTN:  Henry  Cooper 

ATTN:  J.  G.  Lewis 

ATTN:  Tech  Lib 

ATTN:  R.  J.  Port 

Rand  Corporation 

1700  Main  St.,  Santa  Monica,  CA  90406 
ATTN:  C.  C.  Mow 
ATTN:  Tech  Lib 

Science  Applications,  Inc. 

2201  San  Pedro  NE,  Albuquerque,  NM  87110 
ATTN:  J.  L.  Bratton 

Science  Applications,  Inc. 

P.  0.  Box  2351,  La  Jolla,  CA  92038 
ATTN:  Tech  Lib 

Southwest  Research  Institute 

P.  0.  Drawer  28510,  San  Antonio,  TX  78284 

ATTN:  Wilfred  E.  Baker 

ATTN:  A.  B.  Wenzel 

Stanford  Research  Institute 

333  Ravenswood  Ave.,  Menlo  Park,  CA  94025 

ATTN:  Carl  Peterson 

ATTN:  George  R.  Abrahamson 

Systems,  Science  and  Software,  Inc. 

P.  0.  Box  1620,  La  Jolla,  CA  92038 
ATTN:  Donald  R.  Grine 
ATTN:  Tech  Lib 
ATTN:  J.  Sweet 

Terra  Tek,  Inc. 

420  Wakara  Way,  Salt  Lake  City,  UT  87108 
ATTN:  Tech  Lib 
ATTN:  Sidney  Green 


220 


I 


Dept,  of  Defense  Contractors,  Cont'd. 


No.  of 
Co  pi  es 

Tetra  Tech,  Inc. 

630  North  Rosemead  Blvd.,  Pasadena,  CA  91107 
1 ATTN:  Tech  Lib 

TRW  Systems  Group 

One  Space  Park,  Redondo  Beach,  CA  90278 
1 ATTN:  Tech  Info  Center/S-1 930 

1 ATTN:  Peter  K.  Dai  Rl/2170 

1 ATTN:  Norm  Lipner 

TRW  Systems  Group,  San  Bernardino  Operations 
P.  0.  Box  1310,  San  Bernardino,  CA  92402 
1 ATTN:  E.  Y.  Wong  527/712 

TRW  Systems  Group 

Room  712,  Bldg  527,  Norton  AFB,  CA  92409 
1 ATTN:  Gregory  D.  Hulcher 

Universal  Analytics,  Inc. 

7740  W.  Manchester  Blvd.,  Playa  Del  Rey,  CA  90291 
1 ATTN:  E.  I.  Field 

URS  Research  Co 

155  Bovet  Rd.,  San  Mateo,  CA  94402 
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